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Beugung an einer einbackigen Blende endlicher Dicke und der 
Zusammenhang mit der Theorie der Seitenversetzung des total- 
reflektierten Strahles 


j ; Vom Kurt Artmann 
(Mit 1 Abbildung) 


Inhaltsübersicht 


Auf Grund einer an anderer Stelle veröffentlichten Arbeit!) wird die polarisierende 
Wirkung des Blendenschirms auf die Seitenversetzung :des totalreflektierten Strahles 
beim Versuch von Goos und Hänchen?) quantitativ in Rechnung gestellt, was auf eine 
zusätzliche Seitenversetzung herauskommt. Im Anschluß daran ergeben sich einige Be- 
merkungen, zum Teil richtigstellender Art zu einer Arbeit von v. Fragstein?). 


1. Problemstellung 


In einer vorangehenden Arbeit*) hatte der Verfasser die von Goos und Hän- 
chen?) beobachtete Seitenversetzung berechnet, die eine Lichtwelle bei der 
Totalreflexion am dünneren Medium erfährt. Diese Theorie machte nach Kirch- 
hoff die vereinfachende Annahme, daß die Blende, welche die Lichtwelle seitlich 
begrenzt, vollkommen schwarz ist. In der Praxis ist diese Annahme nie erfüllt, 
sondern es müssen auf dem Blendenschirm die Randbedingungen der Maxwell- 
schen Theorie befriedigt werden. 

Eine solche Lösung der Maxwellschen Gleichungen ist von Sommerfeld ) 
für die unendlich dünne einbackige Blende (Halbebene) angegeben worden. 
In konkreten Fällen kann die Sommerfeldsche Theorie daher nur auf einen 
Blendenschirm angewandt werden, dessen Dicke klein gegen die Lichtwellenlänge 
ist. In fast allen praktischen Fällen, speziell beim Versuch von Goos und Hänchen 
ist die Dicke des Blendenschirms jedoch groß gegen die Lichtwellenlänge. 
Deshalb hat der Verfasser an anderer Stelle!) das Beugungsfeld einer Lichtwelle 
hinter einem „dicken“ einbackigen Blendenschirm berechnet. Die wesentlichen 
Ergebnisse dieser Arbeit sollen im folgenden Abschnitt 2 mitgeteilt werden. 


2. Das Beugungsfeld einer Lichtwelle hinter einer einbackigen Blende in der 
Nähe der Schattengrenze 


Die Oberfläche der metallischen einbackigen Blende, welche die Lichtwelle 
beim Versuch von Goos und Hänchen seitlich begrenzt, möge aus einem Zylinder- 


1) K. Artmann, Z. Physik 127, 468 (1950). 

„) F. Goosu. H. Hänchen, Ann. Physik (6) 1, 333 (1947); ferner Ann. Physik (6) 

251 (1949). 

8) C. v. Fragstein, Ann. Physik (6) 4, 271 (1949). 

2 K. Artmann, Ann. Physik (6) 2 ‚87 (1948). ‘ 

Sommerfeld, Math. Ann. 47, 317 (1896). 
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mantel bestehen, dessen Mantellinien senkrecht auf der Zeichenebene stehen und 
dessen Grundkurve die in Abb.1 gezeichnete Gestalt hat. Der Radius a des Halh- 
kreises, der die beiden parallelen Geraden verbindet, soll > Wellenlänge A der 
einfallenden Lichtwelle sein. In praktisch vorkommenden Fällen liegt a zwischen 
10? cm und 102 cm, so daß bei sichtbarem Licht a/A numerisch ~ 10° > 1 ist, 


ebene 
\ 
/rennfläche 


geom. opt 
Schaltengrenze 


geom. opt. 
Schattengrenze 


Abb. 1. Zur Entstehung der Schirmversetzungen 4g, 4 g und der Totalreflexions 


desselben Profils, wie es nach der elementaren Kirchhoff-Huygensschen 
Beugungstheorie zu erwarten ist, an der Schattengrenze um die Strecke 


Ag = 0,20(a/A)'!s -A (1) 


senkrecht zur Einfallsrichtung in Richtung des geometrischen Blendenschattens 
versetzt, falls der magnetische Vektor der einfallenden Welle || Blendenöffnung 
schwingt. Im umgekehrten Polarisationsfall (€ || Blendenöffnung) beträgt die 
Versetzung 
Ag = — 0,39 +A (2) 
\ in umgekehrter Richtung. Die relative Versetzung Ag—Ag bei den beiden Polar 
sationsfällen beträgt demnach 


Arcı = 0,59 (a/A)' + A. (3) 


Wegen a/A>> 1 ist dies eine große Anzahl von Lichtwellenlängen. Mit den zu Be 
ginn dieses Abschnitts angegebenen numerischen Daten ist 


(3) 
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von Goos und Hänchen. Mit kleiner werdender Winkeldifferenz 9,—d, wachsen 
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3. Anwendung von Abschnitt 2 auf den Versuch von Goos und Hänchen 


Wir vergleichen diese Schirmversetzungen 49, Ag mit denjenigen Versetzungen 
Dy, Dg, welche der seitlich begrenzte Lichtstrahl bei der Totalreflexion an 
Se: Trennfläche der beiden Medien 1 und 2 senkrecht zur Fortpflanzungs- 
richtung erfährt (Totalreflexions-Versetzungen). Es seien A, und n, Wellenlänge 
und Brechungsexponent im dichteren Medium (Glas) und A, und n, die entsprechen- 
den Größen im dünneren Medium (Luft); ferner sei 9, der gegen die Normale der 
Trennfläche gemessene Einfallswinkel und #, der Grenzwinkel der Totalreflexion 
(Abb. 1). Dann war nach der Theorie des Verfassers im Falle, daß die magnetische 
Feldstärke $ in der Trennfläche liegt: 

_ 1 Ay 

De= (2) 
während sich im umgekehrten Polarisationsfall (€ in Trennfläche) der Wert 
hy 
Vsin? —sin? 0, 
ergab. In Abb. 1 steht im 1. Polarisationsfall (4g, Dg) der magnetische Vektor 5 
und im 2. Polarisationsfall (Ag, Dg) der elektrische Vektor € der Lichtwelle | - 
Zeichenebene. 

Falls der Einfallswinkel 9, nicht in unmittelbarer Nähe des Grenzwinkels der 
Totalreflexion 9, liegt (9,—d, Z 10°), haben die Totalreflexions-Versetzungen 
Dg, Dg nach (4), (5) die Größenordnung A/z. Nach (3°) ist dann die relative Schirm- 
versetzung Aye) = Ag—Ag rund 10 bis 20mal größer als die Totalreflexions- 
Versetzung Ds, De. In diesem Fall überwiegt die polarisierende Wirkung des 
Blendenschirms gegenüber der Totalreflexions-Seitenversetzung beim Versuch 


(5) 


Dow 


die Totalreflexions-Versetzungen Dg und Dg an, bis sie für d,—d, ~ 1/,° die 
Größenordnung 7-/ der relativen Schirmversetzung Arc = Ag—Ag erreichen. 
Allerdings fällt dieser Einfluß der Blende bei der praktischen Durchführung des 
Versuches nicht ins Gewicht, weil Goos und Hänchen den Lichtstrahl nicht nur 
einmal, sondern siebzigmal und mehr totalreflektieren ließen, während der 
Strahl nur ein einziges Mal die Blende durchsetzt. Auf diese Weise werden die 
Totalreflexions-Versetzungen Dg, Dg mit einem Faktor 70 multipliziert, während 
die Schirmversetzungen Ag, Ag ihren Wert (2), (3) unverändert beibehalten. Aus 
diesem Grunde können bei der Goosschen Versuchsgenauigkeit von 10%, die 
Schirmversetzungen 4g, Ag gegenüber den Totalreflexions-Versetzungen Ag, 
De vernachlässigt werden. 
N 4. Bemerkungen zu einer Arbeit von v. Fragstein ee 
Im Anschluß an die in Anmerkung 4 der $. 209 zitierte Veröffentlichung des 
Verfassers macht Hr. v. Fragstein®) auf eine sehr interessante Arbeit von Schae- _ 
fer und Pich®) aufmerksam, in welcher bereits die Vorgänge bei der Totalreflexion 
einer seitlich begrenzten Lichtwelle behandelt sind. Weder Herrn Prof. Goos q 
noch Frau Dr. Lindberg-Hänchen noch mir ist diese Arbeit bekannt gewesen. = 
Bei Schaefer-Pich wird eine Näherungslösung der Maxwellschen Glei- _ “a 
chungen verwendet, welche der Begrenztheit des Lichtbiindels Rechnung tragt. E 


*) C. Schaefer u. R. Pich, Ann. Physik (5) 80, 245 (1937). 
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pon es wind nicht die Frage behandelt, auf welche Weise diese seitlich begrenzte 
Lichtwelle entsteht. Sollen die Ergebnisse dieser Rechnungen experimentell nach- 
geprüft werden, so muß unbedingt ein Blendenschirm in den Kreis der Betrach- 
tungen einbezogen werden, wie dies hier in den vorangegangenen Abschnitten 2 
und 3 geschehen ist. Es ergibt sich hieraus die Notwendigkeit, die Betrachtungen 
von Schaefer-Pich bezüglich der beiden Polarisationsmöglichkeiten der ein- 
fallenden Lichtwelle, wie im Abschnitt 3 besprochen, zu modifizieren, von denen 
_ Schaefer-Pich nur die eine Polarisationsmöglichkeit behandelt haben (€ || Blen- 

denöffnung). 

Im Anschluß hieran halte ich es für erforderlich, zwei kurze Bemerkungen zur 
_ Arbeit von v. Fragstein zu machen. Einerseits hebt Hr. v. Fragstein hervor, 

daß die Verwendung des Poyntingschen Vektors bei den Schaefer-Pich- 
sehen Betrachtungen besonders anschaulich sei. Bei der elektromagnetischen 
Welle ist eine solche Behandlung zwar möglich. Aber der Effekt der Seitenver- 
setzung des totalreflektierten Strahles an sich bezieht sich ja nicht nur auf elektro- 

magnetische Wellen, sondern auf Wellen beliebiger Art (Schallwellen, Wasser- 
 wellen, Materiewellen). Dann existiert aber der Begriff des Poyntingschen 
_Vektors nicht mehr, während der Begriff des Phasensprungs, den der Ver- 
_ fasser zur Berechnung der Strahlversetzung verwendet hatte”), sich auf Wellen 
beliebiger Art übertragen läßt. 
Die zweite Bemerkung betrifft die Behauptung Hrn. v. Fragsteins, daß die 
Ergebnisse meiner (in Anmerkung 4 der 8. 209 zitierten) Rechnungen bereits voll- 
_ ständig in der Arbeit von Schaefer-Pich enthalten sein sollen. Diese Behaup- 
tung Hrn. v. Fragsteins zwingt mich leider zu einer Richtigstellung, weil sie 
wesentlich über das hinausgeht, was sich aus der Arbeit von Schaefer-Pich ent- 
nehmen läßt. Denn einmal wird — wie Hr. v. Fragstein selbst zugibt — der 


sind bei Hrn. v. Fragstein 4 Seiten Rechnungen erforderlich, um mein Ergebnis 
aus den Schaefer-Pichschen Betrachtungen zu folgern. Und schließlich behandelt 
Hr. v. Fragstein nur den einen Polarisationsfall (€ || Blendenöffnung). Der 
andere Polarisationsfall ist — wie oben bereits bemerkt — nicht einmal von Schae- 
_ fer-Pich behandelt worden, und bei Hrn. v. Fragstein wird für diesen 2. Polari- 
___ gationsfall nur das Ergebnis einer nicht mitgeteilten Rechnung angegeben. 
e 
7) Vgl. Formel (2,8) der in Anm. *) der S. 209 zitierten Arbeit des Verfassers. 
*) Dies geschah erstmalig von Goos u. Hänchen; vgl. Anm. ?), 


Effekt der Seitenversetzung von Schaefer-Pich nicht ausgesprochen ®). Ferner, 
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Berechnung der optischen Konstanten starker PER a 


Inhaltsiibersicht 

Formeln für die Brennweite, die Brennpunktslagen, die Projektivvergrößerung, 
die Farbabweichung und den Öffnungsfehler der starken Magnetlinsen des Über- 
mikroskops als Funktion der Linsenstärke bzw. der Spulenerregung, Beschleuni- 
gungsspannung und Polschuhparameter werden angegeben. Das gegebene Pol- 
schuhfeld wird dazu durch ein ,,Newtonsches Abbildungsfeld“, in denen allein der x 
Begriff der Brennweite sinnvoll definiert werden kann, möglichst gut approxi- 
miert und ein Verfahren zur direkten Berechnung der optischen Konstanten wird 
als entwickelt. 


Zur hohen Vergrößerungen im muß des Ob- 
jekt sehr nahe an das Objektiv herangebracht werden. Es befindet sich daher Ss. 
bereits in Gebieten sehr starker magnetischer Feldstärke und die bekannte Formel?) = 
fiir die Brennweite schwacher Linsen a 

+00 


kann deshalb nicht mehr verwendet werden. Wir haben daher in einer früheren 
Arbeit?) strenge Formeln für die optischen Konstanten für ein typisches magne- 


tisches Abbildungsfeld der Gestalt ’ eras 
B 

B, z = —_! 9 (2 

1+ (2) 


wie es in den magnetischen Objektiven des Ubermikroskops angenähert ver- 
wirklicht ist, angegeben. 

Im Folgenden wollen wir zeigen, wie man die gefundenen Ergebnisse zur Be- 
stimmung der optischen Konstanten empirisch ausgemessener Polschuhfelder ver- 4 
wenden kann. Zur Frage der Feldausmessung vergleiche man die Arbeiten von 2 
0.Klemperer und H. Miller*), L. Marton‘), J. Dosse°) und die in letzter 
Zeit angegebene Methode von M. v. Ments und J. B. Le Poole’). 


1) H. Busch, Arch. Elektrotechn. 18, 553 (1927). 
*) W. Glaser, Z. Physik 117, 285 (1941). ~ 

3) O. Klemperer u. H. Miller, J. sci. Instrum. XVI (1921) 1939. 

*) L. Marton, Physic. Rev. 55, 762 (1939). 

*) J. Dosse, Z. Physik 117, 437 (1941). 

*) M. v. Ments u. J. B. Le Poole, Appl. Sci. Research. Vol. B. 1947. 
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Die Differentialgleichung der achsennahen Bahnen 
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gesetzt worden ist. a bedeutet eine Länge des Wir 
wollen im folgenden dafür — so wie früher — die rechte Halbwertsbreite des Magnet- 
_ feldes B, (z) wählen. B, ist der Maximalwert des Feldes, dessen Abszisse wir als 
 Anfangkpunkt für die Zählung von z wählen wollen. Der dimensionslose Para- 
meter k? bestimmt die „Linsenstärke‘“ d.h. die Stärke des Feldeinflusses auf die 
_ Elektronenbahn, deren Krümmung ihm nach (4) proportional ist. Auch anschau- 
lich ist klar, daß man die „Linsenstärke‘“ durch diesen Parameter zu kennzeichnen 
hat. Denn der Feldeinfluß aut gleich schnelle Elektronen d. h. gleichen U-Wertes 
5 A und bei gleicher Feldbreite a wird um so größer sein je größer B, ist. Ferner bei 
gleichem B, und a um so kleiner je größer U ist, denn einerseits ist der Strahl 
„steifer“ sudsnmeiti sind die Elektronen nur kürzere Zeit dem Felde ausgesetzt. 
Von zwei Feldern mit gleichen .B,-Werten werden gleich schnelle Elektronen von 
dem Felde stärker beeinflußt werden, das die größere Halbwertsbreite besitzt, denn 
der Einwirkung dieses Feldes sind sie länger ausgesetzt. Aus dimensionellen 
Gründen folgt daher, daß man die ,,Linsenstirke“ (bis auf den belanglosen Faktor 
1/,) durch A? in (5) zu definieren hat. Bei sehr hohen Elektronengeschwindigkeiten 
(z.B. U >100kV), muß man die relativistische Massenveränderlichkeit berück- 
N s sichtigen. Man erhält so für (5) 


e Bia? 
8m, U U(1 +525 a vu) 
Jede Lösung von (4) wird k* als Parameter enthalten: 
r=r(x, kt). 


Es sei (7) gerade jene Elektronenbahn, welche achsenparallel in der Kin 
es höhe 1 von rechts einfalle. Die Nullstelle xp 


bestimmt dann den dingseitigen Brennpunkt, mit welchem im Ubermikroskop 
f R praktisch die Dinglage zusammenfällt. Die Brennweite ist nach der bekannten 
Konstruktion durch 


Entwickelt man (7) nach dem Parameter k? der Linsenstärke, so erhält man, da 
fir k* = 0,.r = 1 wird’) 


r= 1—k7,(x) + kr,(a)— (10) 


7) Die wechselnden Vorzeichen wurden aus formalen Gründen gewählt, damit sich 
im Folgenden für die Funktionen r, günstigere Vorzeichen ergeben. 


r (2p, k*) = 0 (8) 


ı 0 
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und für die Brechkraft ergibt sich ayant Sead 

+ (ap) — (ap) (11) 


Der nächstliegende Weg einer wirklichen Berechnung der Brechkraft 1/f besteht 
nun darin, daß man nach (10) die Entwicklungskoeffizienten r, in 


oo 
r= = (- 1y 7,(z), (12) 


aus der Differentialgleichung bestimmt. Wir setzen dazu wie in einer friiheren 
Arbeit ®) (12) in (4) ein und erhalten durch Vergleich gleich hoher Potenzen von k* 


ry = B(x) = 0) (13) 
oder 


Damit kann im Prinzip aus r, = 1 und f(x) die Funktion r,, aus dieser 7, usw. be- 
stimmt werden. 
Man findet so z. B. für 7, 


n= n= Sa f pode. 
oo oo co 
Nimmt man an, daß der Brennpunkt bereits außerhalb des Feldes liegt, so daß für 


alle Werte x< xp die Funktion f (x) praktisch bereits Null ist, so kann man für 
(11) auch schreiben 


(15) 


+00 


(16) 


Das erste Glied dieser Reihe stimmt aber gerade mit (1) überein, wenn man hier 
nach (5) k? einführt. Man sieht also, daß (1) das erste Glied in der Entwicklung der 
Linsenbrechkraft nach der Linsenstärke k®? darstellt. Die Formel (1) gilt somit 
nur für schwache Linsen k? < 1. Da man in der Ubermikroskopie mit k*-Werten 
arbeitet, welche zwischen 1 und 2 liegen, kann somit die Formel (1) zur Berech- 
nung der Brennweite nicht verwendet werden. Man muß also auch die höheren 
Potenzen der Linsenstärke k* berücksichtigen und als obere Grenze in den auf- 
tretenden Integralen xp und nicht — oo wählen. 

Dieses Verfahren — das insbesondere von M. v. Ments und J. B. Le Poole 
verwendet wird — erfordert jedoch einen großen Rechenaufwand. Aus (10) muß 
die erste Nullstelle x = xp bestimmt werden. Dies ist, wenn man sich k* vorgibt 
sehr langwierig, weil es die Auswertung und Tabellierung der Funktionen r, (2) 
für beliebige x-Werte verlangt. Aber auch die umgekehrte Aufgabe der Bestim- 
mung der zu einer vorgegebenen Brennpunktslage xp erforderlichen Linsenstärke k?, 
verlangt die Ermittlung der Nullstellen k) der Funktion 


(ap, = 1 — (tp) + (tp) — (ap) + = 9, (17) 
was bei Beriicksichtigung nur der angeschriebenen Glieder bereits einen relativ 


großen Arbeitsaufwand erfordert, auch wenn man in Betracht zieht, daß im all- 
gemeinen die Glieder r, mit steigendem Index » rasch abnehmen. Schließlich ist 


8) W. Glaser, Z. Physik 81, 451 (1933). 
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das ganze Verfahren grundsätzlichen Bedenken ausgesetzt, worauf am Schluß 
dieser Arbeit eingegangen wird. 

Was man für die praktische Übermikroskopie braucht, sind handliche, geschlos- 
sene Formeln, die ohne nennenswerten Rechenaufwand die optischen Konstanten 
als Funktion der Polschuhdaten, der Amperewindungszahl und der Beschleuni- 
gungsspannung zu bestimmen gestatten. Dabei zerlegt man sich die Aufgabe 
zweckmäßigerweise in zwei Teilaufgaben. Man gibt zunächst die optischen Kon- 
stanten z. B. x, und f als Funktion der Linsenstärke k* an. Dieser Parameter geht 
ja allein in die Differentialgleichung der Elektronenbahnen ein. Nachträglich be- 
stimmt man den Zusammenhang zwischen k* und den anderen Daten. Für 


+00 i 
Ba=p Bidz, (18) 


wobei p eine dimensionslose, allein von den Verhältnissen der Polschuhdimensionen 

abhängige Größe ist. Wir nennen sie die „Gestaltsfunktion“®). Sie kann als 

Funktion der Längenverhältnisse durch Messungen von zusammengehörigen 
+ Werten von B,a und des von der Feldkurve eingeschlossenen Flächeninhaltes für 

4 jeden Polschuhtypus bestimmt werden ?°), Auf Grund des Ampereschen Gesetzes 
i) 


ist nun 
+00 


f B,d&=wnI, (19) 


wobei nJ die Zahl der Amperewindungen der Magnetlinse bedeuten. Wenn wir 
dies in (18) einsetzen, erhalten wir entgültig für die Linsenstärke 


wobei 4,, 62,...die Längenverhältnisse sind, welche die Polschuhdimensionen 
kennzeichnen. Wenn wir praktische Einheiten einsetzen, ergibt sich 


0,276 p (2) 


wobei wir zur Abkürzung De “xe 
U, = U (14+ 0,98- 10 U) (22) 


gesetzt haben. Es bedeutet U, die infolge der relativistischen Massenveränder- 
lichkeit an Stelle von U einzusetzende Beschleunigungsspannung. 
Der Feldverlauf 


sei nun entweder auf Gena von Feldmessungen tabellarisch, zeichnerisch oder als 


analytischer Ausdruck gegeben. Um das ganze von — oo bis + oo für z in Betracht 


®) Eine analoge „Gestaltsfunktion“ (,,Spulenformkonstante“) ist auch in der 
Arbeit von E. Ruska und M. Knoll, Z. techn. Phys. 12,‘ Nr. 8, 389 (1933) ein- 
geführt worden. Vgl. hierzu auch M. v. Ments u. J. B. le Poole, l.c. 

1°) Die Gestaltsfunktionen p(5) für bestimmte Polschuhtypen sind von Wolf 
und Müller in der Abt. f. Tidtecnsncstth der Siemens & Halske A. G. in Berlin- 
Siemensstadt bestimmt worden. Diese Ergebnisse sollen bei anderer Gelegenheit 
veröffentlicht werden. (Zusatz bei der Korrektur.) : a 
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kommende Intervall auf ein endliches Intervall zu reduzieren, führen wir durch 
x=octgp (24) 


einen Hilfswinkel g ein. Diese Substitution wird sich auch in anderer Hinsicht als 
sehr zweckmäßig erweisen. Die Konstante o lassen wir zunächst noch offen. Wenn 
wir (24) in die Differentialgleichung (4) einsetzen, ergibt sich 


cosy dr 


k? 0? B(o ctg ¢) 
=0. (25) 
Setzt man — um die erste Ableitung wegzuschaffen, 
(26) 
sin 
so erhält man En 
4,7 t 
vg) + (1+ (27) 
Aus dieser Gleichung erkennen wir zunächst folgendes: fiir die Feldform 
B(o ctgp) = sin! p (28) 


kann (27) unmittelbar integriert werden. Setzt man außerdem o = 1, so besitzt 
sie in diesem Falle die beiden Lösungen 


w(p) = sin yı +k29, u(p)= cos +Kk:o (29) 


Führt man in (28) statt des Hilfswinkels 9 wieder die Koordinate x ein, so ergibt 
sich 


(30) 


und man erkennt, daß das entsprechende Feld mit dem bereits früher behandelten 
Glockenfeld (2) übereinstimmt. Für die in der Einfallshöhe 1 von rechts achsen- 
parallel einfallende Elektronenbahn, erhalten wir aus (29) und (26) 


1 sin 1+ ko 17 
ya —, z=ctgo. 3l 


Die Nullstelle ergibt sich zu 


= 32 
Yr (32) 
und die ENGE ist daher durch 
2 tp = act 33 
gegeben. Für die Brennweite f folgt aus (9), (31) und (32) 


Das Feld (2) mit den optischen Konstanten (33) und (34) stellt den tatsächlichen 

Feldverlauf in Polschuhlinsen insbesondere im Sättigungsgebiet in den meisten 

Fällen mit genügender Genauigkeit dar"). Im allgemeinen aber haben die ver- 
1) E, Ruska, Arch. Elektrochem. 38, 102 (1944). : 


. 
‘ 
ake 4 
[2 
te ‘ 
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FU 
# 
4 
4 
34) 
= 


Annalen der Physik. 6. Folge. Band 7. 1950 


wendeten magnetischen Abbildungsfelder einen steileren Abfall als er durch (30) 
N ist. So — z. B. dem Kreisstromfeld 


| + (3) Bir 


ct = sin’ 36 


so daß im Unendlichen d.h. für gy —0 die Funktion B(cte@) wie g® gegen Null 
geht. Der Abfall empirisch ausgemessener Felder erfolgt gewöhnlich noch rascher, 
wenn man in einem Feldstärkebereich arbeitet, der von der Sättigung der Pol- 
schuhe noch entsprechend entfernt ist). 

Nun ist das Gebiet des stärksten Feldeinflusses auf die Elektronenbahn wohl 
das in der Umgebung des Feldmaximums für z= 0 d.h. fir p= 5 - Durch 
eine entsprechende Wahl der Konstanten o in (24) können wir nun gerade in diesem 
entscheidenden Feldgebiet, das empirische Feld durch ein Glockenfeld vom Typus 
(30) bis auf Glieder vierter Ordnung in x approximieren. Da 2 = 0 den Ort des 
Feldmaximums bedeutet und das Feld als symmetrisch um x = 0 vorausgesetzt 


werde, hat B(x) die Entwicklung 
Es ist also 
2 
B(octgg) =1+ (38) 
wobei durch die Punkte Glieder vierter Ordnung angedeutet sind. Setzt man nun 
o 2 
so wird 
1 2 1 4 ... 40 
Ploctgp) = 1—2ctg?p + ipa t (40) 
Damit erhält man für die Differentialgleichung (27) Bee 
++ (a) 


Wenn wir daher, die durch die Glieder vierter Ordnung bedingte Korrektur bei- 
seite lassen, erhalten wir die beiden Lösungen 


w(p) = sin yı +ho-p, u(y) = cos Yı+ ko (42) 
und die von rechts achsenparallel in der Einfallshéhe 1 einfallende Elektronenbahn 
durch 
sin 


gegeben. Für die optischen Konstanten erhält man daraus in Analogie zu (3) 
und (34) 


wor 


wer 
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—| 


= — sin ———. 
o V 1+ke 
Die Größe o ist dabei allgemein durch (39) aus dem Feldverlauf bestimmt. Ist 6 


die kleinste Stufe der Abszisse, welche man bei der punktweisen Ausmessung des 
Magnetfeldes lings der Feldachse erreicht, so ergibt sich nach (38) 


2 [1 — 
woraus sich o zu 
V1—B() 


bestimmt. Wenn die Funktion f(x) tabellarisch gegeben ist, bedeutet also f (6) 
den Funktionswert für die kleinste Abszisse. 


Auch aus der Krümmung r der Feldkurve ß (x) im Scheitel kann o entnommen 


werden. Wegen ß’ (0) = 0 ist Be 
= = -#’(0)>0, of= (48) 


Zur Bestimmung von o zeichne man sich auf durchsichtigem Papier eine Schar von 
Kreisen, die alle durch den gleichen Punkt gehen und deren Mittelpunkt auf einer 
Geraden @ liegen. Zu den einzelnen Kreisen denke man sich die Größe o = 2/9 
angeschrieben. . Legt man den Punkt P in den Scheitel von f(x) und läßt die 
Gerade G mit der Ordinatenachse zusammenfallen, so hat man jenen o-Wert zu 
wählen, dessen Kreis sich am besten der Feldkurve f(x) im Scheitel anschmiegt. 
Außer den beiden angegebenen optischen Konstanten (45) und (46) spielt noch 
die Projektivvergrößerung eine wichtige Rolle. Bei dieser in der letzten Stufe er- 
folgenden Vergrößerung handelt es sich um keine eigentliche Abbildurg, sondern 
um eine optische Projektion, so wie etwa bei der „Abbildung“ eines Gegenstandes 
mittels einer punktförmigen Lichtquelle im Schattenbild. Da infolge der hohen 
Vergrößerung die Elektronenstrahlen praktisch achsenparallel in die Ebene des 
Zwischenbildes eintreten, haben wir die von links achsenparallel in das Projektiv 
eintretenden Elektronenbahnen zu betrachten und ihren Achsenabstand in der 
Ebene des Bildschirmes zu bestimmen. Ist g die Einfallshöhe dieses Strahles d. h. 
die ,,GegenstandsgréBe“ auf dem Zwischenbildschirm, so lautet die gesuchte 
Elektronenbahn 
sin 


Man sieht unmittelbar, daß für x = — ood.h.p = die Funktion r den Wert g 
annimmt und man kann sich auch leicht überzeugen, daß hier die Ableitung nach x 
verschwindet, die Bahn also tatsächlich achsenparallel ist. Bedeutet Jp die Ent- 
femung des Endschirmes von der Projektivmitte also —_— die Abszisse 
4,=apctgg, des Leuchtschirmes, so gilt ae 


(49) 


4 
rch (30) und 
(35) 
at 
Ex 
Fr ‘ 
is 
‘ 
(43) 
zu (3) 
(44) 
4 
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Für die Vp=G/g(G =r(@,)) erhält man aus (49) 
g y1+ + sin 9, 
oder wenn man bedenkt, daß !p >apalso @, praktisch gleich Null gesetzt werden 
kann und sin 9, aus (50) ersetzt (valine 
lp sin Vi+ ko PAR 


Vp= 2 
P apo (52) 


Hierbei bedeutet ap die Halbwertsbreite des Projektivfeldes und kp und o beziehen 
sich — wie alle in diesem Absatz verwendeten Größen — gleichfalls auf das Pro- 
jektivfeld. 

Der Radius des auf die Dingebene -bezogenen Offnungsfeblerscheibchens wird 
durch 


(51) 


65 = (58) 
gegeben, wenn a die Strahlapertur ugd C,, die Öffnungsfehlerkonstante bedeutet. 
Für C; ergibt sich in Analogie zu der in der angeführten Arbeit gefundenen Formel: 
G4 zt k? o 148-3 . 22 1 


Hierbei ist eine hohe Vergrößerung vorausgesetzt, bei welcher das Objekt 

praktisch im dingseitigen Brennpunkt angenommen werden kann. Auch den Aus- 

druck für den Farbfehler können wir auf den betrachteten allgemeineren Fall er- 

weitern. Für das auf die Dingebene bezogene Farbfehlerscheibchen ergab sich 


e 


’ wenn AU den Spielraum der Voltgeschwindigkeit bedeutet. Für die Farbfehler- 
konstante erhalten wir 


Ro 1 
V1+ 


In Tabelle 1 haben wir die numerischen Werte der optischen Konstanten in Ab- 
hängigkeit von den maßgebenden Parametern zusammengestellt. Dabei haben 
wir der Vollständigkeit halber auch die dingseitigeHauptpunkts-Abszisse angeführt. 
Wenn wir aus (48) für o einsetzen, erhalten wir z. B. für die Brennweite den voll 
ständigen Ausdruck i 

Wir wollen davon absehen die entsprechenden Ausdriicke auch fiir die anderen 
optischen Größen explizit anzuführen. 

Wenn das empirische Feld z. B. durch einen analytischen Ausdruck der Ge- 
stalt 1%) 
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Tabelle 1 
k? o? Zp/b | b/f za/b C,/b 
0,2 —3,563 0,270 0,138 3,251 15,024 
0,4 —1,890 0,467 0,249 1,734 3,308 
0,6 1,298 . 0,612 0,841 1,246 1,571 ne 
0,8 —0,970 0,718 0,422 1,012 0,999 
1,0 —0,760 0,769 0,495 0,876 0,735 
1,2 —0,609 0,854 0,562 0,792 0,591 
1,4 —0,492 0,897 0,625 0,734 0,502 
1,6 —0,396 0,929 0,679 0,694 0,442 
1,8 —0,316 0,953 0,732 0,664 0,401 
2,0 —0,248 0,970 0,782 0,641 0,370 
2,2 —0,187 0,982 0,830 _ 0,625 0,347 
24 —0,134 0,991 0,875 0,612 0,329 
26 —0,085 0,996 0,918 0,602 0,312 
2,8 —0,041 0,999 0,959 0,598 0,305 iG 
3.0 0,000 1,000 1,000 0,589 0,294 
ko zr/a b/f | zu/b C,/a O,/a 
3,2 0,038 0,999 1,039 0,584 0,287 
3,4 0,073 0,997 1,075 0,581 0,286 
3,6 0,106 0,994 1,113 0,579 0,275 
3,8 0,137 0,990 1,147 0,578 
4,0 0,167 0,986 1,182 0,577 
4,2 0,195 0,981 1,214 0,577 
4,4 0,222 0,976 1,247 0,578 
4,6 0,249 0,970 1,278 0,579 
4,8 0,273 0,964 1,309 0,580 
5,0 0,297 0,958 1,340 0,581 
6,0 0,404 0,927 1,481 0,592 
7,0 0,495 0,896 1,613 0,605 
8,0 0,577 0,866 1,734 0,621 
9,0 0,652 0,837 1,843 0,637 
10 0,719 0,811 1,954 0,653 


Die optischen Konstanten der magnetischen Abbildungsfelder als Funktion der 
charakteristischen Linge b=a0 = | 2B, (0) [B!’(0) bedeutet die zweite Ableitung 


von B,(z) nach z für z= 0] und des Parameters k? o? = - Cp und OC; beziehen 
sich auf hohe Vergrößerung. : 
dargestellt werden kann, erhält man die entsprechenden optischen Konstanten n 


folgendermaßen. Führen wir zunächst die Halbwertsbreiten a gemäß 


ein, so wird nach (5) ER =. 


B} (a _ (60) 
+ - 
und hieraus folgt 


B(2) = 


q 
(49) 
(51) 
werden 
= 4 
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EN 
Vn (se _ 1) 
Die beiden optischen Konstanten vn und (45) erhalten damit die Gestalt 
ctg (63) 
3 
a\2#-1 
ı 


wir 

(2% - 1) 
RR eaten ergeben sich für die anderen Größen. Diese Erweiterungen 
unserer Formeln auf die steiler abfallenden Feldtypen (58) sind im wesentlichen 
bereits von K. SIEGBAHN im Jahre 1943 angegeben worden"). Für den Grenz- 
fall u — co wird 


1-24 (65) 


wobei die Punkte Glieder höherer Ordnung in 1/u andeuten. (60) erhält die Ge- 
stalt 


x) =. lim lim - zu, (66) 
Fiir das Abbildungsfeld vom Typus einer Gausschen Glockenkurve 


B,(r) = Be 


ergibt sich somit für o nach (66) 
1 1 
ij — = 


und die o optischen Konstanten (44), (45) und (52) erhalten dis Gestalt 


(69) 


| » Vin 2 y! +5 


und 
l 


Die 


wal 


gilt 


wo 


Di 


be: 


so daß sich nach (39) fiir o ergi Die 
we 
der 
seht 
Fall 
wir 
geb 
gele 
erfa 
und 
= 
| 
: 
= 
80 
als 
Wi 
Pr fw 
| ™) 
ur P / 
Viti 
“4 14) K. Siegbahn, Arkiv f. Mat. Astr. och Fysik 80, 1 (1944). oe) tee u 


alt 


(64) 


erungen 
ntlichen 
| Grenz- 


(65) 


die Ge- 


(66) 


(67) 


(68), 


(69) 


(70) 


W.Glaser: Berechnung der optischen Konstanten starker magnetischer Elektronenlinsen 223 


Die Länge b hängt dabei mit der Halbwertsbreite a durch die Gleichung 
ben (72$ 


Vin 2 
zusammen. 

Wie die Erfahrung zeigt, sind die Werte der optischen Konstanten zumal die 
der Gaussschen Dioptrik gegen kleine Abänderungen der Funktion f(x) nicht 
sehr empfindlich. Die Formeln (44), (45) und (52) werden also in den meisten 
Fällen für die Praxis bei Dimensionierungsfragen ausreichen. Trotzdem wollen 
wir im folgenden noch eine rechnerisch relativ einfache Korrektur daran an- 
geben, welche die Abweichungen der empirischen Feldkurve von dem zugrunde 
gelegten approximierenden Glockenfeld auch in größerer Entfernung vom Scheitel 
erfaßt. Dazu setzen wir in der Differentialgleichung (4) j 


und erhalten 
qa t Blo (74) 
Die Funktion o? k? ß (o £) hat die 
o? k? Bio £) = o® k2 (0) + = (0) #+--- (75) 


wählen wir o so, daß V 
2 


gilt, so können wir schreiben 
02 k2 


B(o 8) = (77) 
wobei die Punkte Glieder vierter Ordnung bedeuten. Nun schreiben wir statt (74) 
der | 1 1 
de +0? k? 1+ ep" = ker; ß (o 9] r. (78). 
Die Funktion auf der echten Seite von En welche wir mit 0? k?x(£) also 
9 
bezeichnen wollen, beginnt mit der Entwicklung Jat [ 
(80) 


so daß in der Nähe von £=0.d. h. im Gebiet des stärksten Feldeinflusses x (£) 
als klein von vierter Ordnung aufgefaßt werden kann. Für große Werte von & 
wird z(£) wiederum Null. Wir werden daher x(£) lediglich als eine Störungs- 
funktion in der Differentialgleichung für die Elektronenbahnen 


= x(é) (81) 


aufzufassen haben 5). Analog zu einem früher behandelten Fall können wir daher 
die die Lösung mittels Variation der Konstanten bestimmen. Zwei unabhängige Lö- 


15) Über die Zurückführung dieser Differentialgleichung auf eine Integralgleichung 
und ihre Lösung für Feldtypen (58) vgl. J. Schlögl, Wiener Ber. 155 237 (1949). 


i 
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“sungen des ungestörten Problems sind ® 


1 sinop _ 1 coswp 7 
wobei 
@ = yı +ko® und €=ctgp (83) 


gesetzt worden ist. w ist dabei die in der Einfallshöhe 1 von rechts achsenparallel 
_ einfallende Elektronenbahn. Die Lösung von (81), welche den gleichen Anfangs- 
bedingungen genügt, lautet, 


sinop 1, coswp 
Sere sin p (84) 
wenn wir ¢(g) und y(y) durch = 


definieren. Die zum Brennpunkt gehörende Größe 9, = gr bestimmt sich nach 
(84) aus der Gleichung 


[1 + e(pr)] sin —Y(pr) Cos = O. (87) 


(88) 


Wenn man daher 


in die Gl. (87) einsetzt und konsequenter Weise Glieder höherer Ordnung in 7 
vernachlässigt, erhält man 


sin (2 + = —y (=) oder 6=2 (.=(2)), (90) 


so daß sich für die Brennpunktskoordinate z7 = ao ctg 


a+y, 


zp = acct 91 
‚ergibt. Für die Brennweite folgt aus 


or 


$ (84) der Ausdruck 
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dargestellt. Die beiden Korrektionsgrößen y, und &, bestimmen sich nach (85) 
und (86) durch 


: _B (o ctg 
sin! & 


und 


sin 2wad«. 


Die Größe o ist durch un gegeben. Die durch y, und ¢, an den optischen 
Konstanten bewirkte Korrektur wird im allgemeinen nur gering sein. Man wird 
die beiden Integrale (95) und (96) graphisch oder En bestimmen, indem 
man sich die Funktion a 

(o ctg x) 
& (97) 
baw. die Integranden von (95) und (96) tabuliert. 

Zweckmäßiger noch als das angewandte Störungsverfahren zur Berechnung 
der korrigierten Werte der optischen Konstanten erscheint uns die folgende Methode. 
Sie ist im wesentlichen mit der Raleighschen Methode zur genäherten Berech- 
nung der Eigenwerte eines Eigenwertproblems gleichbedeutend. Wir schreiben 
dazu die Differentialgleichung (27) in der Gestalt 


vu’ (9) + AF (gy) v(y) = 9, (98) 
wobei wir F gleich Zu 5 
.2 <2 (0 ctg p) 
F (9) =] k? 0? einig (99) ; 


gesetzt haben. Wir verlangen, daß die Lösung v(p) in py = g und g = q ver- 
schwinde 

=9, =90 (100) 
also z= aoctgg, und 2, =acctggp, einander wie Ding und Bild entsprechen 
mögen. Multipliziert man (98) mit v und integriert die erhaltene Gleichung 
zwischen @, und Pr 80 ergibt sich wegen (100) mittels partieller Integration 


(101) 


Hätte man eine strenge Lösung v(p) von (98), für welche die Randbedingungen 
erfüllt sind, so ergäbe sich daraus der Eigenwert A. Setzt man in (101) für v(p) 
eine angenäherte, die Randbedingungen erfüllende Lösung ein, so ergibt sich 
nach Raleigh der Eigenwert A trotzdem mit einer bemerkenswerten Genauig- 
keit und zwar stets als zu groß. Nun ist 


v = sin w(y —q) (102) 
eine wae Lösung von (98) fiir 2 = 1, die den Randbedingungen genügt, 


Wenn w = gewählt wird. Und zwar ist @? als Funktion von k* o* durch (83) 


gegeben, wenn mi man die Glieder vierter Ordnung in F nicht in Betracht zieht. In 
diesem Falle wäre (101) streng gültig. Setzen wir daher für v(y) aus (102) in 
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(101) ein, erhalten wir: Wipe q 
| 
| ‚#2 cos? w (P — Yo) dp 3 
Nun wollen wir insbesondere die Abbildung des unendlich gre auf den 
Brennpunkt betrachten. Wir wählen dazug, = 0 undg, = —. Damit ergibt sich 


Diese Gleichung gibt uns den gesuchten Zusammenhang zwischen k*o* und o 
Durch Einsetzen erhält man 


208 B octg—) . 
sin! 
0 
Wir sehen, daß für das approximierende Glockenfeld, für welches 
B (octg 
| =] (105) 
sin? — 


ist, w? durch (83) gegeben wird. Allgemein können wir schreiben 


9 B (octg 2) 
g(w) =— | ————sin' gdp. (107) 
sint 
0 


Diese Gleichungen können wir nun durch sukzessive Approximation lösen. Für 
' das approximierende Glockenfeld ist g(w) = 1. Wir setzen daher als erste Lösung 


(108) 
in die rechte Seite von (106) und erhalten : 


=1+ Bo?) (108) 


auf diese Weise kénnte man fortfahren. Es wird aber im es bereits der 
erste Schritt genügen. Damit erhält man ' 
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Die optischen Konstanten sind dann durch 
dargestellt, wo o durch (39) bzw. 6 =ao durch b= ) BO) definiert ist. 


Wir wollen hier ausdrücklich bemerken, daß den Formeln (112) eine grund- 
sätzliche Bedeutung zukommt. Denn diese Gestalt haben ganz allgemein die 
Ausdrücke für Brennpunktslage und Brechkraft in den elektrisch-magnetischen 
Abbildungsfeldern in denen, die gewöhnliche Linsengleichung für beliebige Linsen- 
stärke gilt. In diesen Feldern — die wir ,, Newtonsche Abbildungsfelder“ genannt 
haben — ist es allein sinnvoll von einer Brennweite zu sprechen. Gehört das ver- 
wendete Abbildungsfeld nicht zur besonderen Klasse der Newtonschen Felder, 


so kommt der nach dem Verfahren von Gl. (11) berechneten Größe F keine physi- 


kalische Bedeutung zu, sie ist eine reine Rechengröße. Denn diese Größe erhält 
ja nur dadurch ihren besonderen praktischen Wert, daß sie die Vergrößerung und 
den Zusammenhang zwischen Ding- und Bildort zu bestimmen gestattet, wenn 
die gewöhnliche Linsengleichung gilt. Die Berechnungsmethode durch numerische 
Integration der Elektronenbahn liefert also im physikalischen Sinne die Brenn- 
weite'®) nur ebenso in einem angenähertem Sinne, als sich das zur Abbildung prak- 
tisch verwendete Polschuhfeld durch ein ,,Newtonsches Abbildungsfeld“ appro- 
umieren läßt. Die Formeln (44), (45) und (52) haben wir erhalten, indem wir das 
empirische Feld durch das einfachste Newtonsche Abbildungsfeld möglichst gut 
darstellten. Man kann sich fragen: Läßt sich die Annäherung eines vorgegebenen 
empirischen Feldes durch kompliziertere Newtonsche Abbildungsfelder nicht 
vielleicht noch weiter treiben ? Da die optischen Konstanten aller Newtonscher 
Felder die Gestalt (112) mit zunächst unbestimmtem & haben, kommt praktisch 
die Frage darauf hinaus, den Zusammenhang zwischen w und den das empirische 
Feld kennzeichnenden Daten (z. B. k?, o usw.) für das best-approximierende New- 
tonsche Abbildungsfeld zu bestimmen. Die obige Methode der Eigenwertbestim- 
mung nach Raleigh, die aus der großen Zahl der bekannten Verfahren?’) ver- 
wendet worden ist, gibt vielleicht einen ersten Ansatz dazu, wie die Berechnung von 
® aus dem vorgegebenen Feld ganz allgemein zu erfolgen hat. ae 


Verf. dankt der Vereinigung österr. Industrieller für die Förderung seiner _ 
Arbeiten. 


*) Die Tragweite des Brennpunktsbegriffes und einige hierher gehörigen Fragen 
rden in einer gemeinsam mit O. Bergmann verfaßten Arbeit behandelt, die sich bei der 
leitschrift „ZAMP“ in Druck befindet. 

”) Vgl. hierzu z. B. L. Collatz, Z. angew. Math. Mech. 19, H. 4 u. 5 (1939). Die An- 
wendung der verschiedenen Verfahren der Eigenwertberechnung auf die Bestimmung 
der elektronenoptischen Konstanten soll von einem Mitarbeiter behandelt werden. Vgl. 

auch P. Funk u. W. Glaser, Z. Physik 102, 603 (1936). 


Wien, Institut für angewandte Physik der Technischen Hochschule und Berlin, 
Abt. f. Elektronenoptik der Siemens & Halske A. G. 


3 (Bei der Redaktion eingegangen am 1. Februar 1950.) 
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Es werden z wii 3 nome Methoden zur Berechnung des elektromagnetischen 
Feldes eines translatorisch bewegten, endlich ausgedehnten Körpers gegeben, 
Dann wird gezeigt, wie man vom Standpunkt der speziellen Relativitätstheorie 
aus diese 3 neuen und die 2 bekannten Feldberechnungsmethoden anwenden muß 
bei mehreren, mit verschiedenen Translationsgeschwindigkeiten bewegten Einzel- 
körpern, bei rotierender Materie und bei beliebig bewegten Körpern. Zum Schluß 
werden in 2 Anhängen die erforderlichen Rechnungen und Ergänzungen gegeben. 

Statt der sonst meist benutzten speziellen Lorentztransformation wird 
durchgehend die allgemeine Lorentztransformation verwendet. Ebenso werden 
statt der sonst üblichen speziellen, e und w enthaltenden Materialgleichungen in 
der ganzen Arbeit die allgemeinen Materialgleichungen mit ® und M benutzt. 
Diese beiden Verallgemeinerungen bewirken eine wesentliche Vereinfachung der 
Darstellung und der Rechnungen. 


er I. Einleitung 

Zur Berechnung des elektromagnetischen Feldes bewegter Körper benutzt 
man bisher entweder die Transformationsformeln für die Feldgrößen E, B, H und 
®D (1. Methode) oder die Transformationsgleichungen für die Potentiale g und 
(2. Methode). Formeln zur Berechnung der retardierten Potentiale aus den Raum- 
und Flächenladungsdichten 9 undo, den Raum- und Flächenstromdichten j und jg, 
sowie den elektrischen und magnetischen Polarisationen ® und M (3. Methode) sind 
bisher für in beliebiger Richtung bewegte endliche (begrenzte) Körper noch nicht ange- 
geben worden. Dasselbe gilt für die Ermittlung der momentanen Potentiale aus, 0, 
i, in, B, M und € (4. Methode) und für die Darstellung der Feldgrößen €, ¥B; 9 
und ® aus ihren Quellen und Wirbeln (5. Methode). 

Der Grund für das Fehlen der 3., 4. und 5. Methode liegt in Folgendem: Es 


treten bei diesen Methoden Raumintegrale auf über die Größen = , = 


bzw. ı06®, rot und rot ®. Bei den Integrationen über die Grenzfläche des 


bewegten Körpers sind entsprechende Ausdrücke =, = und Rot € ba. 


Rot ®, Rot $ und Rot ® zu berücksichtigen. Für diese Größen waren aber bisher 
keine Formeln bekannt. In einer Arbeit!) über ,,Die elektrischen und magnetischen 
Flächenwirbel bei bewegten Körpern“ ist nun Folgendes gezeigt worden: 
Irgendein Vektorfeld % werde ' it einem Körper mitbewegt; dann ist an det 
mit der Geschwindigkeit » bewegte Grenzflächen des Körpers für einen ruhenden 


3 1) T. Schlomka, Ann. Physik (6) 5, 190 (1949/50). J 
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Beobachter ein (lokaler) zeitlicher Feldsprung 45 vorhanden, der gegeben ist durch 
A 
= (9:1 — Be) Una: (1) 


Bei Benutzung dieser Formel und Verwendung der daraus resultierenden Flächen- 
wirbel für die elektromagnetischen Feldgrößen lassen sich nun die 3., 4. und 5. 
Methode zur Berechnung des elektromagnetischen Feldes bewegter, endlich aus- 
gedehnter Körper leicht ableiten. 

Wir geben in Abschnitt II der Systematik wegen, und weil wir die darin vor- 
kommenden Gleichungen später brauchen, zunächst die beiden bereits bekannten 
Methoden 1 und 2 und dann die neuen Methoden 3, 4 und 5. Das Bezugssystem 
des Beobachters, für den das Feld berechnet wird, ist ein im Sinne der speziellen 
Relativitatstheorie ,,berechtigtes“, d.h. ein im Fixsternsystem ruhendes oder dagegen 
in irgend einer Richtung mit irgendeiner konstanten Translationsgeschwindigkeit 
bewegtes Koordinatensystem. 

In Abschnitt II handelt es sich um einen einzigen Körper, der dem gewählten 
berechtigten Bezugssystem gegenüber mit einer konstanten Translationsgeschwin- 
digkeit bewegt ist (,,elektrodynamisches Ein-Körper-Problem‘‘) oder um mehrere 
Einzelkörper, die alle dieselbe Translationsgeschwindigkeit haben. Ab- 
schnitt III enthält Ausführungen zum „elektrodynamischen Mehr-Körper-Problem“ ; 
bei ihm sind mehrere Einzelkörper mit verschiedenen Translationsgeschwin- 
digkeiten vorhanden. Abschnitt IV behandelt als Sonderfall des „elektrodyna- 
mischen Unendlich-viel-Körper-Problems“ den rotierenden Körper, und Ab- 
schnitt V den beliebig bewegten Körper. 

Alle Gleichungen sind im rationalen Vier-Grundgrößen-System geschrieben ; 
die Magnetisierung M ist als $ dimensionsgleiche Größe angesetzt. 


II. Elektrodynamisches Ein-Körper-Problem 

A) Erste Methode 
©, B°, D° und H° seien die Felder im ,,Ruhsystem“, d.h. in dem mit dem be- 
wegten Körper mitbewegten Bezugssystem. Der ruhende Beobachter bewegt sich 
dem Ruhsystem gegenüber mit der beliebig gerichteten Translationsgeschwindig- 
keit—p; seine Felder sind dann gegeben?) durch die Transformationsgleichungen : 


E=k (@ - vx B) + (2) 

+ + 0-H (= 9°) = (5) 


nit k = 


a 
= y 
Körper 
f 
4 
| \ 

Viese Gleichungen sind fürs V akuum zuerst angegepven tv. vans vel A. Brill, a 
Vorlesungen zur Einführung in die Mechanik raumerfüllender Massen, Leipzig 1909, S. 217. iF 
Allgemeine Behandlung der Transformationsgleichungen für beliebige Bewegungsrichtung : = 
K:Tamaki, Mem. College of Science and Engineering, Kyoto Imp. Univ. 8, 103, = 
113, 141 (1911/12) und 5, 235 (1912/13); W. v. Ignatowsky, Arch. d. Math. u. Phys. (3) Paes 
li,1 und 18, 17 (1911); C. Kraft, Bull. Int. de l’Acad. d. Sc. de Cracovie, Cl. sc. math. 4 
En (A)1911, 596 und 1912, 385; J. Frenkel, Lehrb. d. Elektrodynamik, Bd.1, 286/297 4 
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=A 
B) Zweite Methode 


Aus den Ruhsystem-Potentialen @° und [° ermittelt man zunächst die Po. 
tentiale 9 und A für den ruhenden Beobachter nach den Transformationsgleichungen 


E | (6) 


é 
€ und B= ihe (8; 9) 
+ 


Die Berechnung von ® und § erfolgt (vgl. Pr. A) nach 
und (10; 11) 
0 
dabei sind ® und M aus PP und M° zu ermitteln nach den Gleichungen 


wx B) (ly 

ra C) Dritte Methode 

Aus den Ruhsystem-Werten 9°, o°, j°, jf, PP und M° erhält man die ent- 
sprechenden Werte fiir den ruhenden Beobachter nach den Gleichungen (12; 13) 
sowie nach 


k) (= (14 
| = k vie v if 7 17 
+2) und o=k(or+ 2h). 


aus berechnet man dann (vgl. Anhang A) die retardierten Potentiale 


fo — di —Di } —rlc 
= + [© is rl df (18) 


, oF 
| rot 
2 vie +(¥%,— Be) Un. + Rot M} ay (19) 
0 


4ur 


und daraus € und ® nach (8; 9), sowie ® und § nach (10; 11). Die Raumintegrale 
sind im allgemeinen über den ganzen unendlichen Raum und die Flächenintegrak 
über die Grenzflächen des bewegten Körpers zu erstrecken. Dabei ist jedoch zu 
beachten, daß der Körper für den ruhenden Beobachter nicht dieselbe Gestalt hat, 
wie im Ruhsystem, sondern in der Bewegungsrichtung die Lorentzkontraktion 


l = 1°/k erfahren hat. af 
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D) Vierte Methode 


Unter Benutzung von (12. . . 17) berechnet man (vgl. Anhang A) die Momentan- 
potentiale 


— div D 
(20) 


1 Pr i+3 rot M+ 
mom = € 
j Rot M — € In 
+ {<= — ¥.) 1,2 + & (&, 2) ® 12 (21) 


und daraus € und ® nach (8; 9), sowie D und § nach (10; 11). Auch hier sind die 
Flächenintegrale über die Grenzflächen des bewegten, durch die Lorentzkon- 
traktion deformierten Körpers zu erstrecken; das Raumintegral in (20) ist im all- 


gemeinen, das Raumintegral in (21) in jedem Falle (wegen des = t) iiber den 


ganzen unendlichen Raum zu nehmen. 


E) Fünfte Methode 


Die zur Quellen- und Wirbeldarstellung von €, 8 und 9 srforderlichen Raum- 
und Grenzflächen-Rotoren und -Divergenzen können der in Abschnitt I zitierten 
Arbeit von Schlomka entnommen werden. Die Quellen von ® sind durch (30) 
und die entsprechende Flächengleichung gegeben; seine Wirbel ergeben sich aus 
(10, 29 und 1) zu: 


rot D = rot ¥—e, — 2 


bzw. Rot D= Rot $ — 5 (Bı — Una: (22) 


Man hat also (bei hinreichend starkem Verschwinden der Feldvektoren im Unend- 


lichen) die Darstellungen: [(23) 

+rotmM D, —D,z) on Rot M 

B= u, -rot dr + [ato a) (24) 


4nr 


D=—grad( 
5 (k= 
rot $ = Rot B—e — vn 


(26) 


div M Div M 1+ In + (Di un, 


Auch hier ist bei den Integrationen die Lorentzkontraktion des ER Körpers 


zu beachten; die 2 und 2 enthaltenden 23 ae sind in jedem Falle 


über den ganzen unendlichen Raum zu erstrecken. 
Aste. 
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F) Vergleich der fünf Methoden 


Man überzeugt sich an Hand einfacher Beispiele leicht davon, daß bei exak- 


ten Berechnungen die 1. Methode allen andern Methoden überlegen ist, 
Rechnungen nach der 2. Methode sind zwar im allgemeinen streng durchführbar; 
doch empfiehlt es sich oft, bei gegebenem 9° und W° erst © und ®° zu berechnen 
und dann mit der 1. Methode weiter zu arbeiten. Exakte Berechnungen mit den 
_ retardierten Potentialen (3. Methode) sind zwar bei bewegten punktförmigen 
Feldträgern (Punktladungen, Dipolen) möglich; bei endlich ausgedehnten Feld- 
trägern sind aber die Auswertungen von (18; 19) der Lorentzkontraktion und 
Retardierung wegen praktisch undurchführbar®). Die 4. Methode ist zur Be- 
rechnung von € prinzipiell nicht verwendbar: man braucht zur Auswertung von 


(21) nicht nur, wie bei ruhenden Körpern, das 3 im ganzen Raum, sondern außer- 


dem noch den &-Sprung &,—€, an der gesamten Oberfläche des bewegten Körpers. 


MN ae Auch zur strengen Berechnung der andern Feldgrößen ist die 4. Methode, selbst 


bei gegebenem €, der Lorentzkontraktion und der über den unendlichen Raum 
zu erstreckenden Integration wegen praktisch nicht brauchbar. Aus demselben 
Grunde versagt die 5. Methode bei exakten Berechnungen. 


Die 4. und 5. Methode können übrigens, da beide Momentan-Darstellungen sind, 
in einander umgerechnet werden. ®, ist unmittelbar gleich ®,. €, läßt sich in €, 
überführen bei Benutzung der in Anhang B bewiesenen Gleichung 


Verwendung derselben Gleichung onion der Quellen-Wirbel-Darstellung von ® 
 liefert D, > D;. Schließlich führt Einsetzen der Quellen-Wirbel-Darstellung von 


Min H, zu H;. 

Bei stationären Zuständen entfällt die Retardierung und werden alle 
_ Zeitableitungen sowie vn. gleich Null. Es liefern dann die 3., 4. und 5. Methode 
(bei Umrechnung von D; in ®, und von 9; in $,) dieselben Feldberechungsvor- 
schriften, die bei v<c denen der 1. und 2. Methode in praktischer Beziehung 


we durchaus gleichwertig sind. 


Zum Schluß dieses Abschnittes sei noch eine prinzipielle Bemerkung zur 4. und 
5. Methode gemacht. Vom rein mathematischen Standpunkt aus sind beide Methoden 
sicher zulässig. Die 4. Methode basiert aber auf der Gl. (44): div Amom = 9, 
und diese Gleichung ist nicht Lorentz-invariant; es ist also vom Standpunkt 
der speziellen Relativitätstheorie aus unlogisch, mit der 4. Methode zu arbeiten. 
Die 4. Methode kann man nun nach den Ausführungen des vorletzten Absatzes 
in die 5. Methode umrechnen; es ist daher auch die 5. Methode aus physikalischen 
Gründen im allgemeinen abzulehnen. Nur dann, wenn die Lorentz-invariante 


Gl. (41) in (44) übergeht, d.h. für = 0, also für statische und stationäre Zu- 


; 3) Bei i im apie stationären Feldern kann man jedoch die von der Lorentz- 
_kontraktion und von der Retardierung herrührenden Schwierigkeiten dadurch umgehen, 
daß man für Pret und Aret andere gibt, die über das unkontra- 


_ hierte Volumen zu erstrecken sind und die Zeit überhaupt nicht enthalten. Dion so abge- 


änderte 3. Methode läßt sich dann in die 2. Methode umrechnen. 
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stände, wird man die (dann mit der 3. Methode zusammenfallende) 4. Methode 
der Momentan-Potentiale und die 5. Methode der Quellen- und atte 
als physikalisch sinnvoll bezeichnen können. R ; 


III. Elektrodynamisches Mehr-Körper- -Problem 


Im Vorhergehenden war zunächst nur das Vorhandensein eines einzigen be- 
wegten Körpers vorausgesetzt. Sind mehrere, mit derselben Translations- 
geschwindigkeit bewegte Einzelkörper vorhanden, so ändert sich jedoch an den 
bisherigen Ausführungen nichts Wesentliches. Die Einzelkörper beeinflussen sich 
zwar gegenseitig: die Ladung des einen Körpers bewirkt Ladungsverteilungen auf 
den andern; die Ströme eines Körpers magnetisieren durch ihr Magnetfeld die 
andern Körper usw. Aber alle diese Einzelkörper haben gegeneinander keine Ge- 
schwindigkeiten, beeinflussen sich also nicht noch zusätzlich infolge ihrer Relativ- 
geschwindigkeiten. In den Gleichungen des 
Abschnitts II ist in diesem Falle unter dem Zh 
Feld % (Kurzsymbol für €, 8, D, 9) die 
Gesamtheit der von den einzelnen Körpern | 
erzeugten Einzelfelder zu verstehen; die 
Flächenintegrationen sind über die Gesamt- 
heit der Grenzflächen aller bewegten Körper 
zu erstrecken. Dieser Fall gehört noch zum 
elektrodynamischen Ein-Körper-Problem:es 


handelt sich eben gewissermaßen um einen Fi Y 
bewegten Einzelkörper mit zahlreichen 
Hohlräumen. 

3 


Anders liegt der Fall, wenn die Einzel- 
körper verschiedene Translationsge- 
schwindigkeiten haben (Abb. 1). Es sei 
zunächst nur K, mit seinem 9°, 6°, j°, jf, 
¥° und M° vorhanden; infolge der Geschwindigkeit », erzeuge er im ruhenden 
Beobachtungspunkt A ein nach Abschnitt II zu berechnendes Feld %,. Ebenso 
liefere ein zunächst nur alleine vorhandener Körper K, infolge seiner Geschwindig- 
keit d, ein Feld %, in A. Bewegen sich gleichzeitig K, mit », und K, mit 0, + d,, 
so ist aber in A nicht etwa das Feld % = % + %, vorhanden. Denn einmal be- 
einflussen sich ja K, und K, schon gegenseitig, wenn sie relativ zu einander ruhen 
(siehe oben). Jetzt haben die Körper aber noch Gegengeschwindigkeiten: K, hat 
K, gegenüber die Geschwindigkeit »,—v, und K, gegenüber K, die Geschwindig- 
keit b.—v,. Infolge dieser Gegengeschwindigkeiten treten in K, und K, neue (der 
Abstandsänderung wegen sich zeitlich dauernd ändernde) elektrische und magneti- 
sche Felder auf, die in K, und K,neue (zeitlich veränderliche) Ladungs- und Strom- 
verteilungen, sowie neue (zeitlich veränderliche) elektrische Polarisationen und 
Magnetisierungen hervorrufen, die %, in %* und %, in %* umändern. Man kann 
also wohl sagen: In A ist ein Feld %* = %* + %* vorhanden, wo %* alleine von 
dem mit d, bewegten K, und %* alleine von dem mit vd, bewegten K, herrührt; 
aber berechnen kann man diese Felder %* und %* auch in den allereinfachsten 
Fällen kaum. Vernachlässigen darf man die von den Gegengeschwindigkeiten 
in K, und K, erzeugten Felder im allgemeinen nicht; sie sind von nema, GréBen- 


Abb. 1. Elektrodynamisches Drei- 


Körper-Problem 
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Noch schlimmer ist die Sachlage beim Vorhandensein von 3 mit verschiedenen 
Geschwindigkeiten bewegten Körpern (Abb. 1), bei denen z.B. K, die Gegen- 
geschwindigkeit vd, — dv, dem Körper K, gegenüber und gleichzeitig die Gegenge- 


_ schwindigkeit »,—v; dem Körper K, gegenüber hat. 


Abschließend läßt sich sagen: Das elektrodynamische Mehr-Körper-Problem 
ist im allgemeinen wegen des Einflusses der Gegengeschwindigkeiten, auf den in 
der Literatur bisher überhaupt noch nicht hingewiesen wurde, praktisch unlösbar. 


IV. Der rotierende Körper 


Bei einem rotieremden Körper liegt ein ,,elektrodynamisches Unendlich-viel- 


Körper-Problem“ vor: Ein bei A’ gelegenes Volumenelement dr’ (Abb. 2) hat die 


Lineargeschwindigkeit = mx t’; irgendein anderes, bei A” gelegenes Volumen- 


element dr’’ hat die Lineargeschwindigkeit », = xr’; es hat also dr’’ dem dy’ 
gegenüber die ,,Rotations-Gegengeschwindigkeit“ 


DRot. Geg. = Ve — 0, =wx(”—r)=mxd. (27) 


Alle Volumenelemente, die auf einer durch A’’ gehenden Parallelen zu tv liegen, 
haben dieselbe Rotations-Gegengeschwindigkeit. Es gibt aber oo? andere Parallelen 
zu tv, die durch andere Volumenelemente des 
rotierenden Körpers gehen. Demnach sind oo 
Rotations-Gegengeschwindigkeiten für den 
willkürlich herausgegriffenen Körperpunkt A’ 
vorhanden. Alle diese Gegengeschwindigkeiten 
werden nach (27) durch den einfachen Aus- 
druck DRot. Geg. = 10 Xb gegeben ; die Gesamt- 
heit der Rotations-Gegengeschwindigkeiten 
für irgendeinen mitrotierenden Kérperpunkt 4’ 
ist demnach so, als ob der ganze Körper um 
diesen Punkt A’ mit tw rotiert. 

Diese Gesetzmäßigkeit der Rota- 
tions-Gegengeschwindigkeiten ermög- 
licht es, das bei rotierenden Körpern vor- 
liegende Unendlich-viel-Körper-Problem in 
praktisch vorkommenden Fällen zu lösen. 
So konnte z.B. in der Arbeit) ,,Relativi- 
tätstheorie und Unipolarinduktion“ das elek- 
trische Feld, das von einer homogen magng 
tisierten, mit tv || Mt rotierenden Kugel für 
einen ruhenden Beobachter erzeugt wird, 
Abb. 2. Rotations-Gegengeschwin- unter Beriicksichtigung der Rotations-Gegen- 


digkeit geg =%2—01 = geschwindigkeiten nach der speziellen Rela- 


tıvitätstheorie berechnet werden. Vernach- 


ae lässigt wurden dabei die Lorentzkontraktion und der Einfluß der Beschleu- 


nigungen. Die erste Vernachlässigung ist in allen praktisch vorkommenden Fällen 
wegen 9 <c zulässig und muß bei Anwendung der speziellen Relativitäts- 


theorie auf rotierende Körper zur Vermeidung der sonst auftretenden bekannten 
_ Schwierigkeiten gemacht werden. Der Einfluß der Zentrifugalkräfte auf die mit- 


4) T. Schlomka, Ann. Physik (6) 5, 51 (1949). 
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rotierenden freien Elektronen kann nach den Rechnungen von Swann®) (On the 
magnetic and electric fields which spontaneously arise in a rotating conducting 
sphere) im allgemeinen ebenfalls vernachlässigt werden. 

Auf die Bedeutung der Rotations-Gegengeschwindigkeiten wurde schon 1933 
in einer Arbeit ®) über „Gravitation und Erdmagnetismus“ hingewiesen. Der Ge- 
danke der Rotations-Gegengeschwindigkeit ist aber bereits in den Edlundschen 
„Untersuchungen über Unipolarinduktion, atmosphärische Elektrizität und 
Polarlicht“ von 18787) und in seiner „Theorie der Unipolarinduktion“ von 18877) 
enthalten. 


V. Der beliebig bewegte Körper 
Sein Bewegungszustand ist gekennzeichnet durch eine momentane Trans- 
jationsgeschwindigkeit v, und die momentane Winkelgeschwindigkeit tv. Der 
Einfluß der Beschleunigungen tv x [m r’], —* 2 und = xt ist entweder zu vernach- 


lässigen oder gesondert zu berechnen. Die A v, und tv herrührenden Felder % 
können nach Abschnitt II und IV ermittelt werden. Es sind also wegen des tw 
auch beim beliebig bewegten Körper die Rotations- a 
digkeiten zu berücksichtigen. 

Zum Schluß sei nochmals ausdrücklich darauf hingewiesen, daß allen vor- 
stehenden Ausführungen ein im Sinne der speziellen Relativitätstheorie „be- 
rechtigtes“ Beobachtungssystem zu Grunde liegt. Auf die bei Anwendung 
eines rotierenden Bezugssystems vorliegenden besonderen Verhältnisse 
kommen wir in einer andern Arbeit zurück. 


Anhang A. Elektrodynamische Potentiale bei bewegten Körpern 


Nach Min kowski gilt auch bei bewegten Körpern das Maxwell- -Gleichungs- 
system 


roth =j+ = - div) = 0; (28; 30) 


rt —&, div 8 =0. 9; 31) 


Als Materialgleichungen verwenden wir nicht die meist in der Relativitäts- 
theorie benutzten Gleichungen . 


= © = & © und = = Ho 9°. (32; 33) 


Denn einmal sind sie viel zu speziell: sie setzen im Ruhsystem isotropes Ma- 
terial voraus und schlieBen dort permanente Magnetisierung und Elektrisierung 


5) W. F. G. Swann, Terrestrial Magnetism 22, 149 (1917). 

6) T. Schlomka, Gerlands Beitr. z. Geophysik 38, 357/406 (1933); vgl. dort S. 385, 
Abschnitt F 7. 

7) E. Edlund, Kongl. Svenska Vetenskaps-Akad. Handlingar 16, Nr. 1 (1878) und 

22, Nr. 5 (1887); Philos. Mag. (5) 6, 289 (1878) und (5) 24, 401 (1887); Ann. chim. 
physique (5) 16, 49 (1879) und (6) 11, 145 (1887). — Edlund denkt sich nicht, wie 
wir oben (vgl. Abb. 2), den mitrotierenden Aufpunkt A’ ruhend und die Gesamtheit der 

ellpunkte um ihn mit tw rotierend, sondern faßt umgekehrt jeden einzelnen Quell- 
punkt A” als ruhend und den Aufpenkt a als um u einzelnen — A” mit 
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aus. Zum andern führen sie zu unzweckmäßigen Materialgleichungen für den 
ruhenden Beobachter: Einsetzen der zu (2 . . . 5) reziproken Gleichungen in (32; 33) 
liefert 


v 


1—k/» v 1—k/»v 


Skalarmultiplikation beider Gleichungen mit » gibt = EN he 
Einsetzen dieser Werte in (34; 35) liefert die bekannten, recht unhandlichen 


Materialgleichungen 
D+ und B-2xE =" (9—vxD). (36; 37) 


7 Viel zweckmäßiger ist es, von den ganz allgemein gültigen Ruhsysten- 
Materialgleichungen 


De = Bund + Me) (38; 39) 
auszugehen. Setzt man in diesen Gleichungen fiir die darin vorkommenden 6. 


_ VektorgréBen die zu (2...5; 12; 13) reziproken Transformationsgleichungen ein, 
so erhalt man mit 


(40 
und den 
= + P und 
und 4 


mit den einzigen Lösungen © = 0 und R = 0. Es gelten also (bei ruhenden und 
auch bei in beliebiger Richtung translatorisch bewegten Körpern !) für den ruhenden 
Beobachter die handlichen “Materialgleichungen 


wud B=w(P+M). (10; 11) 
Aus (28... 31; 40) folgen in bekannter Weise mit der Nebenbe- 
dingung 
10 
div Get = — (41) 
die Wellengleichungen fiir und Aret: 
a N — = — (2) 


Amom = 0 
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die Poisson-Gleichungen für Pmom und Amom en 


2 Apmom = —div®) (45) 


AG nom = — po (i + rot M + (46) 


Die Lésungen von (42; 43) sind unter den bekannten Voraussetzungen iiber das 
Verschwinden der dabei auftretenden Flächenintegrale im Unendlichen in (18; 19) | 
gegeben, die Lösungen von (45; 46) entsprechend in (20; 21). Man sieht, daß zur 
Aufstellung der über die Grenzfläche des bewegten Körpers zu erstreckenden 
Flächenintegrale in (19) und (21) die Kenntnis von (1) erforderlich ist: 


AB _ b = (€, — 
* = (Bi — Be) On, , ( 2) 


Alle Gleichungen dieser Arbeit gehen natürlich mit » = 0 in die für ruhende 
Körper gültigen Beziehungen über. Die aus (18; 19) und (20; 21) mit vn. = 0 
folgenden Potentialpaare widerlegen also die übliche®) Behauptung, es sei das 
Integrationsverfahren der elektrodynamischen (retardierten oder momentanen) 
Potentiale auf den Fall eines einheitlichen ruhenden Mediums beschränkt. 
Wenn man nicht von vornherein einschränkende Bedingungen macht, also nicht 
die Spezialgleichungen (32; 33), sondern die allgemein gültigen Materialgleichungen 
(38; 39) benutzt, erhält man die für beliebiges ® und beliebiges M geltenden | 
Gleichungen (18... 21). 

Es sei noch ausdrücklich darauf hingewiesen, daß in (18; 19) bei den Retar- 
dierungen c die Vakuum-Lichtgeschwindigkeit bedeutet. Im Gegensatz dazu tritt 
bei den sonstigen Darstellungen *) der retardierten Potentiale für ruhende einheit- 
liche Materie bei der Retardierung die Zeit t — r/c* auf, wo c* = (eu)* ist. Dieser 
scheinbare Widerspruch klärt sich folgendermaßen auf: Die Potentiale sind nur 
Rechengrößen zur Ermittlung von € und ® nach (8; 9). Bei homogener ruhender 
Materie, die den ganzen Raum erfüllt und den speziellen Materialgleichungen 
(32; 33) gehorcht, kann man eben drei verschiedene Potentialpaare auf- 
stellen: Erstens ein Potentialpaar (Pmom ; Amom), bestimmt durch (20; 21), jedoch 
ohne die dortigen Flächenintegrale ; zweitens das beı Hund angegebene Potential- 
paar (fer; Are), bei dem die ,,Wirkungskugel“ sich mit der Geschwindigkeit 
c* = (eu)-+ kontrahiert; und drittens ein Potentialpaar (Pret; Aret), gegeben 
durch (18; 19), jedoch ohne die dortigen Flächenintegrale, bei dem die Wirkungs- 
kugel sich mit der Vakuum-Lichtgeschwindigkeit c = (&14)* zusammenzieht. 
Diese 3 Potentialpaare sind zwar voneinander verschieden, geben aber alle dasselbe _ 
€ und B. 

Daß die Gleichungen (42; 43) auch bei bewegten Körpern gelten, ist bereits 
1911 von Ishiwara?) nn worden. Der Übergang von (38; 39) zu (10; 11) 
wird bei ihm aber überbaupt nicht begründet, weder in der obigen Weise, noch 
durch einen Hinweis auf die Lorentz-Invarianz der Vierergleichung 


Fa- M 


8) A. Sommerfeld, Elektrodynamik, Wiesbaden (1948/49), S. 147; J. Fischer, 
Einf. in die klass. Elektrodynamik, Berlin 1936, S. 87/88. 
®) F. Hund, Einf. in die theor. Physik, Bd. 2, Leipzig 1947, S. 259/265. 
10) J, Ishiwara, Proc. Tökyö Math.-Phys. Soc. (2) 6, 72 (1911). i 
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mit Ay = (9: cD), Fx = (cB; ©) und My = (M;—c $). Die Angabe von 
Ishiwara, die Gleichungen (10; 11) entsprächen den im Ruhsystem geltenden 
re Gleichungen (38; 39), genügt ohne nähere Begründung nicht; denn mit demselben 
Recht könnte man den (falschen) Schluß ziehen, bei bewegten Körpern entsprä- 
chen den Gleichungen (32; 33) die Materialgleichungen D=e€ und =n 
_-wahrend doch in Wirklichkeit die davon verschiedenen Gleichungen (36; 37) für 
den ruhenden Beobachter gelten. Im übrigen gibt Ishiwara auch nur die Dif- 
 ferentialgleichungen (42; 43), aber nicht ihre Lösungen (18; 19); dazu fehlte ihm 
an die hierfiir erforderliche Kenntnis von (1). 
Die Integraldarstellungen des elektromagnetischen Feldes, die 1911 von 
_C. Kraft (vgl. das Zitat in Abschnitt II A) gegeben worden sind, befriedigen in keiner 
_ Weise. Es handelt sich bei seinen Gleichungen (40) und (45) um eine Darstellung 
der Feldtensoren Fy, und Hz, durch Integrale über den unendlichen vierdimen- 
. sionalen Raum. Zu ihrer Auswertung für einen beliebigen Weltpunkt (2, Yo, Zp, ty) 
muß die raumzeitliche Verteilung der Materie-Vierergeschwindigkeit w, Vierer- 
___ stromdichte s, elektrischen Leitfähigkeit o, sowie von e und u für alle Zeiten t< 1, 
gegeben sein. Abgesehen von den Schwierigkeiten einer Auswertung der Kraftschen 
x Integrale i in konkreten Fällen ist unter anderem Folgendes gegen sie einzuwenden: 
Zunächst liegen ihnen nicht die allgemein gültigen Materialgleichungen (10; 11), 
> a die speziellen Gln. (36; 37) zu Grunde; ferner setzen sie voraus, daß die 
elektrische Leitfähigkeit o von Null verschieden ist; schließlich kommen in ihnen 
nur Integrationen über vierdimensionale Volumenelemente dd’ = dx -dy -dz-icdt 
vor, aber nicht über (dreidimensionale) Flächenelemente. Es handelt sich bei 
Kraft also nur um kontinuierliche bewegte Materie; im Gegensatz dazu denkt er 
Er, sich aber (S. 617) ,,die Welt mit einer beliebigen Anzahl beliebig bewegter Körper 
Ei: > ER die irgend welche Gestalt und Größe und verschiedene magnetische 


bilität haben können“ und schließt, damit der Vierervektor Hik - ver- 


hi schwindet, „das Vorhandensein von endlichen Raumintervallen und endlichen 
_ Zeitabschnitten mit stetig veränderlicher Permeabilität aus“. Die Einzelkörper 
sollen also an ihren Grenzflächen Sprungstellen der Permeabilität haben; es 
_ hätten also unbedingt die Sprungflächen der bewegten Einzelkörper berücksichtigt 
werden müssen! Von einer vierdimensionalen Sprungflächengleichung, die unserer 

_ dreidimensionalen Gl. (1) entspräche, ist aber bei Kraft nirgends die Rede; seine 
_ Integraldarstellungen gelten also, wenn sie überhaupt zulässig sind, nicht bei 
endlich ausgedehnten bewegten Körpern. 


Anhang B. Ableitung einer bei bewegter Unstetigkeitsfläche geltenden Formel 
Im Innern eines durch die ruhende Hülle H begrenzten Raumes V befinde sich 
ein Körper K mit seiner Grenzfläche Gr. K sei Träger eines in K und außerhalb 
von K stetigen Vektorfeldes %. Ist dieses auch an der Grenzfläche Gr stetig, so 
gilt bekanntlich die Gleichung 


Wenn g de an Sali unstetig ist, so zerlegt man den Raum V in den von K einge- 
nommenen Raum V, und den Restraum V,, wendet dann (47) an auf V, (mit seiner 


Begrenzung Gr) und auf V, (mit seinen Begrenzungen Gr und H), addiert die so 
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%— 
r 


ist. Man er- 


erhaltenen 2 Gleichungen und berücksichtig, daß n, , x 
hält so die beim Vorhandensein einer Unstetigkeitsfläche Gr geltende 


Formel: 
rot, Bar = + ay + f xd. (48) 
v v Gr H 


Partielle Differentiation nach ¢ lieiert: 


rot, (a, [ea + + (49) 
H 


Bei der Umformung der linken Seite von (49) ist für den Fall, daß K eine Trans- 
lationsgeschwindigkeit » hat, Folgendes zu beachten: K schleppt als Träger von % 
dieses Feld % mit; es treten daher an ruhenden Raumpunkten (lokale) zeitliche 
Änderungen von $ auf. Diese sind im allgemeinen stetig veränderliche Funktionen 
der Zeit; nur an den ruhenden Raumpunkten, über die die Grenzfläche gerade 
hinwegstreicht, springt % plötzlich von einem Wert %, auf den davon verschiedenen 
Wert %;. Nach dem Anhang der in Abschnitt I zitierten Arbeit von Schlomka 
über „Die elektrischen und magnetischen Flächenwirbel bei bewegten Körpern“ 

ist bei Integration über den ruhenden Raum V, der die — Unstetigkeits- 


flache Gr enthalt: 
t 
a | 5 dt 


Setzt man (50) mit (1) in die linke Seite von (49) ein und beachtet, daB das letzte 
Integral der rechten Seite von (49) bei hinreichend starker Abnahme von u nach 
außen verschwindet, so erhält man die bei bewegter Unstetigkeitsfläche 


Gr geltende Formel: 2 


Gr 


[ee] 


(50) 


Herrn Dr. H. Epheser danke ich für einige klärende Besprechungen des 
Themas. 


Hannover, Wilhelm-Busch-Straße 7. ed 


(Bei der Redaktion eingegangen am 31. Januar 1950.) 
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Uber eine statistische Teilfrage zum Ordnungsproblem 
der 2-dimensionalen bindren Mischkristalle 


Die in der Arbeit „Zum Ordnungsproblem in Mischkristallen“:) erwähnte 
Frage nach der Häufigkeit durchlaufender Linien in einem quadratischen Geraden- 
netz, dessen Kreuzungspunkte sich auflösen, wird statistisch behandelt und für 
‘2 numerische Fälle beantwortet. 


Bei einer Untersuchung über die Herstellung der Ordnung in einem binären 
Mischkristall durch Platzwechsel benachbarter Atome, über die in einer früheren 
Arbeit?) berichtet wurde, war die Frage aufgetaucht, ob in einem ebenen quadrati- 
schen Gitter, das mit 2 Atomsorten im Verhältnis 1:1 besetzt ist, bei völlig unge- 
ordnetem Zustande Grenzlinien zwischen ,,richtig“ und „falsch“ besetzten Ge- 
bieten mit merklicher Häufigkeit vorkommen, die das ganze Gitter von einem 
Rande bis zum gegenüberliegenden durchlaufen. Denn solche Grenzlinien werden 
‘durch die nachgewiesene Verkürzungstendenz bei sinkender Temperatur nur aus 

geglättet und quergestellt aber nicht beseitigt, so daß ihr 
| Vorhandensein der Einstellung einer das ganze Gitter um- 
fassenden Ordnung entgegensteht. In der erwähnten Arbeit 
wird die Frage auf das folgende rein statistische Problem zurück- 


geführt. Gegeben sei ein ebenes Netz aus zwei sich gegenseitig 
| . kreuzenden Linienscharen, in dem sich die Linien einer Schar 
nirgends kreuzen, und die Zahl der Kreuzungspunkte auf 
} jeder Linie groß sein soll. Nun sollen sich alle Kreuzungs- 


punkte in der Weise auflösen, daß die vier in ihm zusammen- 
laufenden Linien in der Art zweier Hyperbeläste auseinander- 
treten, wasin jedem Kreuzungspunkt auf zwei Weisen möglich 
ist (Abb. 1). Gefragt wird, mit welcher Wahrscheinlichkeit 
Anordnungen entstehen, die keine durchlaufenden Linien von 
lösungsmöglich- einem Rande des Netzes zum gegenüberliegenden enthalten. 
keiten einer Es ist keine Beschränkung der Allgemeinheit, wenn man die 
Kreuzung Linien jeder Schar gerade, parallel und senkrecht zu denen der 
andern Schar annimmt. Die Zahl der Kreuzungspunkte auf 
jeder Linie sei M. Wird jedem Kreuzungspunkt eine bestimmte der beiden 
Auflösungsarten beigelegt, so heiße eine solche Anordnung eine Komplexion. 
Es gibt also 2% Komplexionen. In jeder Komplexion kommen nur Linien 
_ vor, die entweder geschlossen sind, oder an zwei Randpunkten endigen. Gefragt ist 


1) I. Runge, Ann. Physik (6) 7, 129 (1950). 
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also nach der Anzahl der Komplexionen, in denen keine Linien vorkommen, 
deren Randpunkte auf gegenüberliegenden Rändern des Gitters liegen. 

Wir suchen zunächst die Zahl der möglichen durchlaufenden Linien von einem 
Randpunkt aus. Führen wir ein Koordinatensystem x, y parallel den Gitter- 
geraden etwa mit dem Nullpunkt im linken unteren Eckpunkt ein, und nennen wir 
die Strecke zwischen benachbarten Kreuzungspunkten einen Schritt, so bestehen 
alle Linien aus abwechselnden Schritten ||y und x, die wir als Vor- und Rück- 
wärts, bzw. Seitwärtsschritte bezeichnen wollen, indem wir den Anfang der Linien 
auf der x-Achse etwa an der Stelle 2 = a annehmen (Abb. 2). 

Eine durchlaufende Linie hat mindestens M Vorwärtsschritte, und wenn sie 
M + m Vorwärtsschritte hat, so muß sie m Rückwärtsschritte haben, im ganzen 
sind es also. M + 2 m Schritte ||y. Zwischen je zweien davon liegt ein Seitwärts- 
schritt, also im ganzen auch M + 2m, wenn wir mit einem Schritt ||y beginnen 
und mit einem ||z aufhören. Eigentlich müßte dieser letzte Schritt als längs des 


Randes verlaufend weggelassen werden. Da aber M groß ER ee 
sein soll, ist der Fehler, der durch die Hinzunahme dieses Y r u a 7 
letzten Schrittes begangen wird, sicher klein, und die Formeln mane = 


werden auf diese Weise iibersichtlicher. Die Anzahl der 
durchlaufenden Linien von 2(M + 2m) Schritten wäre y 


demnach das Produkt aus der Zahl F(M, m) der Möglich- i eS 
keiten, mit M+ m Vorwärts- und m Rückwärtsschritten 
die Entfernung M zu erreichen, ohne die Gebiete <O und De R oe 
>M zu berühren, und der Zahl Z(M, m) der Möglichkeiten, u ar 
mit M + 2m Seitwärtsschritten, die =O sein können, ; 

von a aus eine Entfernung zwischen O und M zuerreichen, 9 a M x 


ebenfalls ohne die Gebiete <O und > M zu berühren. Abb. 2. L 
. 2. Lage des 

Die so zu berechnende Zahl ist allerdings wesentlich zu Netzes im Koordi- 
groß, denn sie zählt alle.durchlaufenden Linien mit, die sich natensystem 
selbst wiedertreffen und dann Teile ihres Weges mehr- 
fach durchlaufen. Wegen der Entstehungsart der Linien durch Auflösung der 
Kreuzungspunkte sind solche Linien für unser Problem nicht zulässig, denn 
wenn eine Linie einen bereits berührten Punkt wiedertrifft, so steht ihr als 
nächster Schritt nur noch die Vierte der in dem Punkt zusammenlaufenden Strecken 
frei, während an einem noch unberührten Punkt stets zwei Möglichkeiten der 
Linienfortsetzung vorhanden sind. Es ist nun aber nicht notwendig, die so redu- 
tierte Zahl der Linien selbst zu berechnen, denn in unserm Falle kommt es nicht 
auf diese an, sondern auf die Zahl der Komplexionen, in denen solche Linien vor- 
kommen. Eine durchlaufende Linie, die sich selbst nicht wieder trifft, berührt 
AM + 2m) Punkte. Für diese Punkte ist daher durch die Wahl der Linie die 
Auflösungsrichtung der Kreuzung jeweils in einem bestimmten Sinne festgelegt. 
Bei den übrigen M?—2(M + 2m) Punkten sind noch je zwei Auflösungsrichtungen 
möglich. Die Zahl der Komplexionen, die zu einer sich nicht wiedertreffenden durch- 
lufenden Linie mit gegebenem m.gehören, ist daher 24*-2(M+2m), Betrachten 
wir nun eine durchlaufende Linie, die sich in einem Punkte selbst wieder trifft. 
Von den sämtlichen Fortsetzungsmöglichkeiten die von diesem Punkte an noch 
vorhanden sind, fallen dann die Hälfte fort, weil an diesem Punkt statt zweien 
tur noch ein Schritt frei ist. Aber da die übrigbleibenden möglichen Linien einen 
Punkt weniger berühren, so ist ihre zugehörige Komplexionenzahl zweimal so 
&oß, folglich ändert sich die Komplexionenzahl durch das Fortfallen der Hälfte 
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der Linien nicht. Dasselbe gilt für jeden weiteren Punkt, in dem sich die Linie 
selbst wiedertrifft. Obwohl wieder die Hälfte der dann noch möglichen Linien- 
fortsetzungen ausscheiden, kommen dafür den beizubehaltenden Linien je die 
doppelte Komplexionenzahl zu, weil ein weiterer Punkt von der Linie nicht be- 
rührt wird, also noch zwei Auflösungsmöglichkeiten besitzt. Demnach gilt für 
die Zahl D(M, m, a) der Komplexionen, die von x = a ausgehende durchlaufende 
Linien mit 2(M + 2m) Schritten enthalten, F(.M, m) - E(M, m) - 22°-2(M+2m) 
und fragen wir nach dem Bruchteil der Komplexionen, der solche Linien enthält, 
so ist durch 24° zu teilen, man hat also 

D(M,m,a) = F(M, m) E(M, m) 


92(M+2 m 


Es sind nun die Funktionen F(M, m) und E(M, m) auszurechnen. 


Bereehnung des Faktors F(M, m) 


Unter F(M,m) war die Zahl der Wege zu verstehen, mit M + 2m Vorwärts- 
oder Rückwärtsschritten die Entfernung M zu erreichen, ohne dabei die Punkte 
<0 oder > M zu berühren. Fragen wir zunächst nach der betreffenden Zahl 
ohne diese beiden Bedingungen und nennen etwa f(M,n) die Zahl der Wege, mit n 
Schritten die Entfernung M zu erreichen, so ist f(M, n) gleich der Summe der 
Anzahlen der Wege, mit »—1 Schritten bis M + 1 oder bis M—1 zu gelangen. 
Es gilt also 


Nun gilt für das Pascalsche Dreieck der Binomialkoeffizienten die Rekursions- 
formel 


Drücken wir hierin die Ordnungszahl k der Glieder einer Zeile durch den Ab- 
stand x der betreffenden Stelle von der Mitte der Zeile aus, setzen also in der n. Zeile 


2k—n=x und demnach k= 
so wird die Rekursionsformel > 
n—1 n—1 n 
a 


n—1 a—1}] ' In-ı z+1] In x 


Sie ist dann mit der obigen Formel für f(M, ) identisch, wenn man 2 = M und 


= ae 


setzt. Da nun die erste Zeile des Pascalschen Dreiecks 1, 
1 offenbar die Zahlen /(1,-+ 1) der Wege darstellt, um mit n—= 1 Schritt nach 
M=-+1und M=-—1 zu gelangen, so gilt die Übereinstimmung von f mit den 
Binomialkoeffizienten für das ganze Dreieck. Man erkennt,dies auf anschauliche 
Weise, wenn man etwa jede Ziffer des Dreiecks durch Striche mit den beiden 
schräg darunter stehenden Ziffern verbindet, und sich überzeugt, daß jede Ziffer 
die Zahl der Wege angibt, die längs solcher Striche von der Spitze des Dreiecks 
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bis zu der betreffenden Stelle führen. (Striche nur abwärts zu durchlaufen) (Abb. 3). 
Setzt man noch für die Zahl der Schritte n = M + 2m ein, so ist also 


Wir stellen jetzt die Forderung, daß keine Stellen mit negativem Abstand, 
also links von der Mittellinie des Dreiecks berührt werden dürfen. Hierdurch 
fällt nicht nur die ganze linke Hälfte des Dreiecks (Abb. 3) fort, sondern auch 
noch alle die Wege, die 
von den Ziffern der ersten 
Vertikalen links von der 
Mittellinie ausgehen m 
wieder ins positive zu- — 
rücklaufen. Da nun von 
jeder Stelle aus genau 
gleichviel Wege nach 
beiden Seiten ausgehen, 
so erhalten wir das, was 
an einer beliebigen Stelle 
xin Abzug zu bringen 7 2 7 1 


ist, indem wir in der- Abb. 3. Pascalsches Dreieck der Binomialkoeffizienten 
selben Zeile die zu x in x 
bezug auf die Vertikale <= — 1 symmetrische Ziffer aufsuchen, also z. B. 


für die in der Abbildung durch einen Kreis bezeichneten ‚Stelle die durch ein 
Quadrat bezeichnete Ziffer, die wegen des symmetrischen Baus jeder Zeile mit 
der Ziffer bei x + 2 identisch ist. Denkt man sich die linke Seite des Dreiecks um 
die Vertikale x = — 1 nach rechts umgeklappt und von jeder Ziffer des unteren 
Blattes die auf dem oberen Blatt an derselben Stelle gelegene Ziffer abgezogen, 
(im Beispiel also 15—6 = 9), so hat man jeweils die Zahl der Wege, die das 
Negative nicht berühren. Für die Entfernung M ist also diese Zahl = 


(x +m+ (x zr) m . 


Sollen nun die Wege auch die Punkte rechts von 2 = M nicht berühren, so 
haben wir in derselben Weise das zu subtrahieren, was zu der betreffenden Stelle 
symmetrisch mit Bezug auf die Vertikale x = M + 1 liegt, oder es ist die rechte 
Seite des Dreiecks um diese Gerade nach links umzuklappen und wieder die auf 
dem oberen Blatt stehende Ziffer von der jeweils auf dem unteren Blatt stehenden 
abzuziehen. Da nun aber wegen der ersten Faltung auch von den Ziffern rechts 
von x = M+ 1 etwas abzuziehen ist, so hat man, wenn man sie zunächst ganz 
abzieht, diesen Betrag wieder zuzufiigen. Fiir die ne M ergibt sich 


demnach als Zahl der erlaubten Wege ru 

M+2m M+2m)\ | 


Dies gilt allerdings nur solange m < M + 2 ist; für größere m kommen noch 
weitere Faltungen hinzu, so daß man schließlich für beliebig große Schrittzahl die 
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bei der für jedes endliche m die Glieder von einem gewissen k an verschwinden, 
Bei der Berechnung zeigt sich, daß die Glieder mit k > 0 praktisch schon nicht 
mehr ins Gewicht fallen. Für das Glied k = 0 empfiehlt sich die Umrechnung 


Berechnung des Faktors E(M,m) 


Dieser Faktor zählt die Möglichkeiten, mit M + 2m Seitwärtsschritten von 
einem Anfangspunkt aus Punkte zwischen den Entfernungen —a und M—am 
erreichen, ohne diese beiden Grenzpunkte zu überschreiten (Abb. 2). Ohne die 
letztere Bedingung wären also die Zahlen der (M+ 2m) ten Zeile des Pascal- 


schen Dreiecks von (x. sac “ui 


tm-3 M+m—z 
der Wege über —a nach links und M —a nach rechts hinausgehen, so sind wieder an 
jeder Stelle die Werte der Spiegelpunkte an den beiden Geraden 2 = —a—1 
und « = M-—a-+ 1 abzuziehen und die der Spiegelpunkte dieser Spiegelpunkte 
wieder hinzuzufügen usw. mit abwechselnden Vorzeichen, bis die Spiegelpunkte 
M + 2m 


.) zu summieren. Soll aber keiner 


die Entfernung M + 2 m überschreiten. Schreiben wir zur Abkürzung( yy 


= h(k), so ergibt sich 


M-a 2M+2-a 3M+4-a 


2 


ony 3 2 


Da die Gesamtsumme der (M + 2m) - ten Zeile des Pascalschen Dreiecks 24+?" 
ist, können wir E auch erhalten, indem wir von 2%+2™ das doppelte aller negativen 
Glieder der obigen Entwicklung abziehen, also 


E(M,m,a) = 2 a) * (k) + 

( (*) 4 ) 


Weil M eine große Zahl ist, werden die Zahlen h(k) mit guter Näherung durch 
äquidistante Ordinaten einer Gausschen Fehlerkurve dargestellt und die Summen 
in der Entwicklung können durch Integrale dieser Kurve über entsprechende Inter- 
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valle ersetzt werden. So ergibt sich mit der Abkürzung M+2m=n e 
E(M ‚m,a) 


a> 


—2* dx verstanden wird. 


wenn D(x) das Integral — 


Damit haben wir eine Formel ie die Zahl der Komplexionen, in denen vom 
Punkt a ausgehende durchlaufende Linien der Schrittzahl 2(M + 2m) enthalten 
sind, nämlich F(M, m): E(M,m,a) + 24°-2(M+2m) und für den Bruchteil, 
den diese unter allen 2% Komplexionen bilden, daher 


Summiert man dies über alle m von 0 bis oo, so hat man den Bruchteil K (a) aller 
Komplexionen, die überhaupt von a ausgehende durchlaufende Linien enthalten. __ 


Es ist also 1— K (a) der Bruchteil, in dem keine von a ausgehenden durchlaufenden 
Linien vorkommen. Unter diesen können aber solche mit durchlaufenden Linien, 
die von anderen Punkten ausgehen, enthalten sein. Falls man annehmen darf, 
daß sie unter diesen in der gleichen relativen Häufigkeit vorkommen wie unter 
sämtlichen Komplexionen, so wird daher der Bruchteil ganz ohne durchlaufende 
Linien, der mit W, bezeichnet sei, das Produkt aller (1 — K(a)) sein, also 


Ww,= II (-Ka)= II (-Kü)). 
a=0'--M a=0:--M/2 


Hier sind nur die durchlaufenden Linien berücksichtigt, die auf den Rändern 
zendigen. Für die Rander||y ergibt sich im Fall eines Quadrats noch einmal der- 
selbe Ausdruck. Also wird die Wahrscheinlichkeit, durchlaufende Linien weder 
vom x- noch vom y-Rande aus zu erhalten W = W;. 


Numerische Auswertung D120,mı 


Um eine Vorstellung von = yi 
der Größe dieser Bruchteile u a-M/2 
bekommen, wurden die Formeln tot 
für die Zahlwerte M = 20 und 
50 angenähert ausgewertet. Es 
ergaben sich für die F, E, D 
die in Tabelle 1 angegebenen 
Werte. 

Aus diesen Werten wurden 
die in Abb. 4 und 5 wieder- arr 
gegebenen Kurven als Funk- 
tionen von m gezeichnet und 0 Ww 60 80 100. 120 40 m 


graphisch integriert. Die So g€- = Abb. 4. D(M,m,a) als Funktion von m ‘für 
fundenen Werte sind in 2 M = 20 M d M 
= a=>-— und — 


a=M 
03} 


: 
4 
245 A 
. 
u w 
on 
% 
os 
> 
4 
: 
# 
4 


Annalen der Physik. 6. Folge. Band 7. 1950 


Tabelle 1 M : 
M = 20 aus 


E(M,m,) 


F(M, m) 
“QM+2m 


M M 
QM +2m IM +2m oa, m, am) D{ 


10 | 0.344 - 103 0,7730 | 0,5537 0,266 -10-% 0,1905-10-8 
20 | 0,896-10-* | 0,6070| 0,4299 0,544 -10-3 0,386 -10-3 
30 | 1,116-10 | 0,4750 | 0,3357, 0,530 -10-3 0,375 -10-3 
40 | 1,099-10 | 0,3706 | 0,2621 » | 0,406 -10-% 0,288 -10-3 
60 | 0,832-10= | 0,2270 | 0,1606 | 0,189 -10-% 0,1335-10-3 
100 | 0,393-10- | 0,0836 | 0,0619 0,0329-10- 0,0243 -10-3 Lin 
M = 50 
0,0025-10- | <0,97 | <0,77 | 0,4489 | 0,001 104 | 
0,036-10- | 0,9172| 0,6899 | 0,3894 | 0,033 10-4) 0,025 -10-4| 0,014- pr 
0,364 10-4 | 0,7722| 0,5532) 0,3020 | 0,281-10- 0,201 +10 0,110- die 
0,708 | 0,6364; 0,4512| 0,2446 | 0,449- 10-4) 0,318 -10-4 0,172- pol 
0,836 0,5232| 0,3702) 0,2003 | 0,438-10- 0,310 -10-4 0,167- fal 
0,825-10-* | 00,4390} 0,3033| 0,1640 | 0,362-10-4| 0,250 -10-4 0,135- 
0,720-10-* | 0,8526| 0,2488| 0,1347 | 0,254-10-4| 0,179 - 10-4 0,097- du 
0,560-10- | 0,2360| 0,1670| 0,0908 | 0,132-10-) 0,0935-10-| 0,051: de 
0,238 | 0,1090| 0,0779) 0,0441 | 0,031-10-) 0,022 -10-) 0,0125-10" | H; 
pu 
Tabelle 2 un 
et‘ 
i SO! 
 M=20 K(=) =2,6-10% | M=650 = 1,21- 10-8 
M M 
I—)= . 2 —)= 
K(7) 1,8 -10- K(>) 0,87 -10- 
K (=) = 0,45 - 10+ 
W 
= K(0) =0,0 -10- ge 
Zur Berechnung von W, kann man, da K(a) klein gegen 1 ist, setzen 


D (50m) M=50 Aus den berechneten Werten 
it “re von K(a) läßt sich graphisch 
a=Mi2 der Verlauf von K(a) ange- 
Summe wiederum durch das 

M 
Integral [ K(a)-daersetzen; 

6 

att ergeben sich für M = 20 
und 50 Kurven, die sich 
nahezu nur um einen kon- 
m m stanten Faktorunterscheiden 


_ Abb. 5. D(M,m,a) als Funktion von m für M = 50, (Abb. 6). Wertet man diese 
MoM M 


Flächeninhalte graphisch 
2 Ss) ee so findet man fir 


M M 
2 = - K(a)) = K(a 
: 
+ 


),001 + 104 
),014 104 
),110 104 
),172 
),167: 104 
),135 
097 ‘10° 
,051 107 
‚0125-107 
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M = 20 und 50 die in der 2. Spalte von Tabelle 3 angegebenen Näherungswerte, 
aus denen W, und W bestimmt wurden. 


Tabelle 3 
M nW, | We W=Wi| 1-W 
20 | 0,337 | 0,71 | 0,608 | 0,497 
50 | —0,385 | 0,68 | 0,46 0,54 


Wir finden also, daß der Bruchteil der Komplexionen ganz ohne durchlaufende 
Linien in beiden Fällen nahe bei 50%, liegt, also ist auch die Häufigkeit von An- 
ordnungen mit durchlaufenden Linien etwa ya, 

50%. Die Zahlen erheben nicht den An- f 
spruch, exakte Auswertungen zu sein, da 7% 
die Kurven mit Hilfe ausgiebiger Inter- 2 
polationen gezeichnet wurden; aber jeden- , 
falls zeigt sich, daß durchlaufende Grenzen 
durchaus in einem erheblichen Bruchteil Me Ma MR na 
der Fälle vorkommen; daß ferner ihre pp, 6, Komplexionenbruchteil mit 
Häufigkeit von der Zahl M der Kreuzungs- durchlaufenden Linien von «=a 
punkte pro Gitterseite nur wenig abhängt, aus als Funktion von a 

und daß kaum anzunehmen ist, daß sie 

etwa für sehr große M verschwinden sollte, da sie von M= 20 bis 50 
sogar zuzunehmen scheint. 


Zusammenfassung 

In einem quadratischen Geradennetz sollen sich alle M? Kreuzungspunkte so 
auflösen, daß je 2 sich nicht schneidende Kurvenstücke entstehen. Es wird die 
Wahrscheinlichkeit berechnet, daß eine Anordnung ohne von einem Rand zum 
gegenüberliegenden Rand verlaufende Linien entsteht. Für M = 20 und 50 er- 
geben sich dafür Wahrscheinlichkeiten von ungefähr 50%. 
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Der Einfluß e einer Blende i in Rohren 
ons das Feld einer einfallenden elektromagnetischen Welle 
Von Ernst Ruch 
Inhaltsiibersicht 


In der EHER Arbeit wird das Problem der Beugung und Reflexion einer 
elektromagnetischen Rohrwelle an einer senkrecht zur Rohrwandung angebrachten 
ebenen Blende beliebiger Gestalt theoretisch untersucht. Die Art des Ansatzes 
erfordert die Lösung eines Gleichungssystems mit unendlich vielen Unbekannten, 
Es gelingt, die reziproke Matrix zu konstruieren und damit die Lösung explizit 
hinzuschreiben. Eine näherungsweise numerische Berechnung wird durchgeführt 
für den Fall eines Rohres mit rechteckigem Querschnitt und schlitzförmiger Blende 
bei einfallender H-Welle. Das Problem der Auflösung eines unendlichen Gleichungs- 
systems — wie es sich aus einer gemischten Randwertaufgabe ergibt — wird völlig 
ohne Spezialisierung auf das Rohrwellenproblem diskutiert. Lae 


Einleitung 


Wird in ein ideal leitendes, in Richtung seiner Achse unendlich ausgedehntes, 
zylindrisches Rohr beliebigen Querschnitts aus dem Unendlichen ein elektro- 
_ magnetisches Wechselfeld bestimmter Frequenz eingestrahlt, so bildet sich im 
Innern ein Zustand aus, der sich als Überlagerung von Eigenzuständen auffassen 
läßt. Man unterscheidet diese nach E-, H- und L-Wellen, je nachdem die Kom- 
ponente des magnetischen oder des elektrischen Feldes in Richtung der Rohr- 
achse oder beide verschwinden. Ledinegg!) hat gezeigt, daß unter Voraussetzung 
eines Querschnitts, der im Sinne der Mathematik ein einfach zusammenhängender 
Bereich ist, nur E- und H-Wellen existieren können und aus diesen beiden sich 
jeder mögliche Zustand des Feldes linear kombinieren läßt. Die Aufstellung der 
_ Eigenfunktionen scheint uns besonders durchsichtig mit Hilfe der Debyeschen 
Potentiale und wird daher in dieser Form zu Beginn unserer Abhandlung gebracht. 
Eine senkrecht zur Zylinderwandung angebrachte ebene Blende ruft Beugungs- 
yi 5 und Reflexionserscheinungen hervor, die einer Behandlung schwer zugänglich 


angenommen wird. 

‘ Unter dieser Voraussetzung ist die Berechnung des Feldes mehrfach in Angriff 
genommen worden, jedoch mit der Beschränkung auf sehr kleine Öffnungen oder 
_ sehr kleine Schirme, eine Spezialisierung, die als nullte Näherung die Lösung des 
geschlossenen bzw. völlig offenen Rohres nahe legt. 


1) E. Ledinegg, Ann. Physik (5) 41, 537 (1942). 
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Wir lassen im folgenden ein beliebiges Öffnungsverhältnis zu und überlagern 
der einfallenden Welle ein Wechselfeld, das als Linearkombination von Eigen- 
lösungen mit unbestimmten Koeffizienten angesetzt wird. Das Gesamtfeld muß 
den Randbedingungen auf der Blendenebene genügen, und ihre Formulierung 
im Sinne des kleinsten Fehlerquadrats führt auf ein Gleichungssystem mit unend- 
lich vielen Unbekannten). 

Bei der Auflösung begnügen wir uns zunächst mit einem Fall, bei dem keine 
Umwandlung auftritt, d.h. das modifizierte Feld durch E- oder H-Wellen allein 
beschrieben werden kann. Von den Koeffizienten dieses Gleichungssystems lassen 
sich Eigenschaften nachweisen, die die Konstruktion einer reziproken Matrix in 
Gestalt eines unendlichen Produktes erlauben. Obgleich die Grenzmatrix dieses 
Matrizenproduktes nicht explizit angegeben werden kann, eignet sich ihre Dar- 
stellung zur numerischen Berechnung der Lösung. Eine graphische Veranschau- 
lichung des Feldverlaufs auf der Blendenebene gibt einen Begriff von der Leistungs- 
fähigkeit der Methode. 

Die rein mathematische Auflösung des Systems durch Konstruktion einer in- 
versen Form wird viel allgemeiner durchgeführt, als für unseren Spezialfall nötig 
wäre. Es läßt sich zeigen, daß die Fälle, in denen Umwandlungen auftritt, einer 
Behandlung nach der gleichen Methode zugänglich sind, wenngleich das Ausmaß 
der dafür nötigen Rechenärbeit die Möglichkeit einer vereinfachenden Umformung 
wünschenswert erscheinen läßt. Auch scheinen Probleme außerhalb des Fragen- 
komplexes der Rohrwellen so formulierbar zu sein, daß sie in völlig analoger Weise 
lösbar werden. Wir bringen die Diskussion darüber am Schlusse der Arbeit. 


I. Allgemeine Ableitung der Rohrwellentunktionen 


Ein elektromagnetisches Feld bestimmter Frequenz w wird beschrieben durch 
die Maxwell-Gleichungen in der bekannten Form: 


Die hier angenommene Zeitabhängigkeit ei? soll konsequent beibehalten werden. 
Wir gehen ähnlich wie Debye®) bei der Behandlung der Beugung von Licht 
an einer Kugel von den beiden Ansätzen aus: ne 


E bzw. H 
9=-" wot €=*°" vot 


Daraus folgt durch Einsetzen ih die zweite bzw. erste Mawxell-Gleichung mit 


| 


€=KI/T+ grad® =k IT + grad. 
3 A. Sommerfeld, Vorlesungen iiber theoretische Physik, Bd. VI, Leipzig 1945, 


D ode 


8) P. Debye, Ann. Physik 30, 57 (1909). ae 
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Weiter ergibt sich aus den Maxwell-Gleichungen für den Fall E und H in gleicher 
Weise die Differentialgleichung: 


grad ® — tot rot IT + "®iT=0. 
Ein Koordinatensystem u, v, z sei so gewählt, daß wu, v irgendein Orthogonal- 


system in der Querschnittsebene des Zylinders und z die dazu senkrechte Koor- 
_ dinate parallel zu seinen Mantellinien bedeutet. 


‘ Es empfiehlt sich die Annahme, der Vektor IT habe überall die Richtung der 
positiven z-Achse, also: 


Ii = I, = 37, 
wenn 3 den Einheitsvektor in dieser Richtung bezeichnet. Damit vereinfacht sich 
die Komponentendarstellung der Differentialgleichung wesentlich. Denn wegen: 


gadd=(u 2 rot (V2 ,—u 0) 


tot tot = ( =) U: (7 Uv Lal ar) + (7 
— U und J sind dabei Funktionen von % und 2, die sich aus den speziellen Eigen- 
schaften des u, v-Systems ergeben — gilt: 


dz \U a)” 

‚fayUom , a(V aly) , 


U/V hängt von z nicht ab. Die beiden ersten Gleichungen lassen sich daher durch 

den Ansatz ® = = befriedigen, und man erhält die skalare Gleichung: 4 


_ Mit /\ sei dabei der auf w, v transformierte Laplace-Operator in der Ebene er- 
klart. 

Fiir die elektrischen und magnetischen Feldkomponenten ergibt sich nun im 

einzelnen: 
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Oder in der eleganteren Darstellung des Vektorkalküls und bei Unterscheidung 
zwischen Transversal- und Longitudinalfeld: 


E H 
grad Ge, = “CF vot = —*2" 5 x grad IT 
6, = 
9, 9. kell + . 


An der Rohrwandung miissen die Tangentialkomponenten des elektrischen Feldes 
verschwinden, und daraus folgt: 


E IT=0 H oe 


= = Differentiation in Richtung der Zylindermantelnormalen. 


Wir sehen, die Einführung der beiden Debyeschen Potentiale schließt die 
Trennung in die beiden Wellenformen der Z- und H-Wellen schon in sich. Weiter 
stellen sich die Lösungsfunktionen in ihrer Abhängigkeit von w und v als die Lö- 
sungen der schwingenden Membrane heraus, wie schon Lord Rayleigh*) gezeigt 
hat, und zwar im E-Falle mit festen und im Falle H mit freien Rändern. Selbst- 
verständlich ist die Ableitung nicht an die Orthogonalität von uw und v gebunden. 

Bezeichen wir die Eigenfunktionen der Membrane mit g, bzw. h,, so sind die 
Rohrwellenpotentiale g, und h, gegeben durch: 


E Jr (u, v,2) = 9, (u v) e+ .H hy (u, v, z) =h, (uv) e+ 


Die Eigenwerte der Membrangleichung L; bzw. K; genügen der Beziehung: 
baw. Kij=k-—k. 


Daraus läßt sich schließen, daß /,, k, nur reell oder rein imaginär sein können. Die 
Rohrwellen sind also völlig ungedämpft fortschreitende W. ellen — homogene Wellen 
— oder ein rein abklingender Vorgang — inhomogene Wellen —. 

Die g, und die h, sind jeweils unter sich orthogonal. Sind entartete Lösungen 
dabei, so wollen wir diese im weiteren Verlauf als orthogonalisiert annehmen. 
Außerdem seien sie alle normiert. 


5 
II. Ansatz für das durch eine Blende modifizierte Feld — Er 3a 


einer einfallenden Rohrwelle 


Das Metallrohr soll jetzt an der Stelle z = 0 eine unendlich dünne, ideal leitende 
Blende enthalten. Diese definiert eine Querschnittsebene, die, in ihrer Ausdehnung 
vom Zylindermantel begrenzt, als Blendenfläche angesprochen und durch den 
Buchstaben @ bezeichnet werden soll. Ihre Schnittkurve sei J’, der vom Metall 
erfüllte Teil der Blende heiße Schirm A, daß Gebiet der Öffnung B und der Schirm- 


4) Lord Rayleigh, Philos. Mag. 48, 125. 


Bae 


} 
« 
a 
| 
+ 
x 
u 
3 | 
i 


rand y (Abb. 1). Die Gestalt der Öffnung sei willkürlich und dasselbe gelte für 
den Zylinderquerschnitt. 

E Eine aus z = —oo einfallende E- oder H-Welle vom Typus einer Eigenfunk- 

_ tion erleidet am Schirm Beugung und Reflexion. Wir beschreiben die Veränderung 


_ des ursprünglichen Feldes durch Überlagerung eines Zusatzfeldes. 
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| 
Es sei: 


Cg, bzw. Ch, das Potential der einfallenden Welle mit 
der Amplitude C 


K = g™ + hyp) das Potential für das Zusatzfeld im 


Raum z <0. 
 Abb.1. Gestalt der 
Blende. D = g) + hipy das Zusatzpotential für z > 0. a 

iy oe Die Trennung in zwei Summanden bedeute hier die 


Unterscheidung zwischen einem E- und H-Potential. 


Auf der Blendenfläche z = 0 haben wir einige Rand- und Stetigkeitsbedin- 
gungen zu erfüllen, und zwar: 
1. Verschwinden des elektrischen Transversalfeldes auf dem Schirm: 


rot (h® 3) 


(1a) 
iou 


+ — grad g) + — rot auf A, 


a 
Cz grad 9, 
tot (hy 5) 


c 


[7 

{R) 
+ = grad g(® + 
gradg, 
2. Stetigkeit des elektrischen Transversalfeldes im Loch: 


twp 
c 


rot 3) = gad g(P) + FOR vot (A) 3) auf B. (lb) 


c 


é 
grad + 


3. Stetigkeit des magnetischen Transversalfeldes im Loch: 

F rot (5) + = rot 5) + 2 gradh auf B. (Ie) 
ee Im folgenden wird der Fall einer eingestrahlten Z- bzw. H-Welle dadurch gleich- 
zeitig behandelt, daß in den Formeln die sich unterscheidenden Glieder in der aus 
(la) ersichtlichen Weise übereinander geschrieben werden. 

Denken wir uns g und h nach Eigenfunktionen entwickelt in der Form: 


Er 


G, = S dy” h, 


= Seg, eribz AD) — Sd j, erikez 


alle positiv zu nehmen. Denn für reelle Werte sind damit homogene Wellen ge- 
LEN zeichnet, die von der Blendenfläche weggehen, und im anderen Falle ist die End- 
lichkeit oder genauer das Verschwinden im Unendlichen gesichert. 
Führen wir die Reihendarstellung in die Bedingungen (1) ein, so folgt aus (1a) 
und (1b): 


( 
v 
= 


E.l 


auf 


2. 
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auf A und B. Das bedingt die Relationen: 


wenn wir annehmen, daß alle Summanden in @ zueinander orthogonal sind. 
Dies läßt sich aber mit Hilfe der Greenschen Sätze leicht zeigen. a 


G r G r 
| tot (hy 3) ot (hy. 3) df = | gradi, - gradh, df = ds + Ki Ösu (3) 
6 G r 


J grad 9, tot df = (grad 9, - gradh,) = J tot (9, gradh,) 3 df 


= [9 gradh, -d3 . 


6), bezeichnet das Kronecker-Symbol. Wegen der Eigenschaften von g, 


und h, als Lösungsfunktionen der Membranengleichung fließt daraus die obige 
Asnahme. 


Aus den Beziehungen (2) zusammen mit (1c) folgt das uniinieten des ma- 
gnetischen Transversalfeldes im Loch. 


Wir können jetzt zeigen, daß damit alle zu fordernden Stetigkeitsbedingungen 
für das Feld auf der Blendenfläche erfüllt sind. 


Das Verschwinden des zusätzlichen magnetischen Transversalfeldes im Loch 
zieht das Nullwerden und damit die Stetigkeit der z- Komponente des elektrischen 
Zusatzfeldes innerhalb dieses Gebietes nach sich. Mit dem Verschwinden des ge- 
samten transversalen elektrischen Feldes auf dem Schirm ist Verschwinden und 
Stetigkeit des magnetischen Longitudinalfeldes verbunden. Die Longitudinal- 
komponente des magnetischen Zusatzteldes im Loch ist wegen d(® = d(D) stetig. 
Wir haben Grund anzunehmen, daß das Problem damit eindeutig charakterisiert 
ist, es sei denn, der Schirmrand y läßt Mehrdeutigkeiten zu. 


Es ist wegen der physikalischen Existenz einer solchen Lösung zu erwarten, 
daß Reihenentwicklungen in der benutzten Form sinnvoll sind, d.h. daß sie die 
Funktion überall mit Ausnahme der Schirmkante darstellen. Diese Forderung 
der Darstellbarkeit oder der fast gleichmäßigen Konvergenz solcher Reihen wird 
sicher durch gewisse Stetigkeitseigenschaften an der @hirmkante y gegeben sein. 


Es ist leicht einzusehen, daß der elektrische oder magnetische Kraftfluß durch 
irgendeine Fläche endlich bleiben muß. Das heißt aber Endlichbleiben von Kurven- 
integralen der Form: 


$E-ds 69:48. 


Denken wir uns eine derartige Kurve senkrecht durch die Schirmkante, im 
übrigen willkürlich verlaufend, so folgt daraus insbesondere die Integrierbarkeit 


(la) » 
auf A, 7 

7 
zleich- 
er aus. 
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End- 
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von &, 9, G, 5, nach v über y weg, wenn wir das Koordinatensystem so wählen, 
daß v = const die Schirmkante wiedergibt. 

Da die Integrierbarkeit die Stetigkeit des Integrales in sich schließt, läßt sich 
zeigen, daß die an der Schirmkante tangentiellen Komponenten des elektrischen 


es damit zu begründen, daß wir die Blende als Grenzfall einer metallischen Platte 
endlicher Dicke ansehen können, auf deren Stirnfläche diese Grenzbedingungen 
ja zu fordern wären. Die schlechteste Konvergenz wird also für die auf y senk- 
rechten Komponenten zu erwarten sein. Die Tatsache ihrer Integrierbarkeit 
schließt Konsequenzen in bezug auf das Verhalten ihrer Entwicklungskoeffi- 
- zienten in sich, die wir im folgenden als Konvergenzforderungen an unser Glei- 
chungssystem stellen werden müssen. Die Forderung einer Lösung mit jener 
 Konvergenzeigenschaft wird das Problem eindeutig machen. 

Unsere beiden Postulate, Verschwinden des gesamten transversalen elektri- 
schen Feldes auf A und Verschwinden des magnetischen transversalen Zusatz- 
 feldes auf B, lassen sich im Sinne der kleinsten Abweichung im quadratischen Mittel 
stellen. In Formeln: Es soll 


tot (h, 3) 


4 
und 


gu einem Minimum gemacht werden. = ie 


" Beschränken wir uns zunächst auf endliche Summen, so können diese beiden 
Forderungen, abgesehen von unwahrscheinlichen Spezialfällen, nicht gleichzeitig 
_ einzeln für sich erfüllt werden. Denn denkt man sich die Quadrate der Absolut- 
beträge als Produkt von den konjugiert komplexen Ausdrücken hingeschrieben 
und die Ableitungen nach den konjugiert komplexen Koeffizienten c* und d 
gleich null gesetzt, dann bekommen wir doppelt so viele Gléichungen wie Unbe- 
kannte. Wir multiplizieren daher das zweite Integral mit einer positiven Zahl 
und stellen fiir ein Minimum von dem Ausdruck: 


Man erhält damit endlich viele Gleichungen mit ebenso vielen Unbekannten. 
Die Lösung wird von der Wahl des willkürlich wählbaren Faktors A abhängig sein. 
Wenn aber im Grenzfall n >00 die beiden Integrale für sich verschwinden, und 
das können wir vermuten, da ein vollständiges Orthogonalsystem zugrunde 
N gelegt ist, dann wird die Lösung von A unabhängig sein. A hat die Bedeutung eines 

_ Gewichtsfaktors, da er mit wachsendem Betrag für eine endliche Anzahl von 
Gleichungen das Verschwinden des zweiten Ausdrucks der gleichen Forderung 
an den ersten gegenüber bevorzugt veranlaßt. Wir werden auf die Unabhängigkeit 
von A nach einer age? Diskussion des Gleichungssystems noch später 


y 


und magnetischen Feldes €, und §,, stetig sein müssen. Am leichtesten jedoch ist, 


ul 


di 
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rot (g d, grad df | 
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‚wir nachweisen, daß es den Energiesatz bereits enthält. Zu diesem Zweck ist die 
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Durch die angedeutete Variation nach cy und dy ergeben sich die Gleichungen : 


li l, Cy | grad In grad 9, df d, 7 grad 9, tot (hy 3) df 


Ä 
+4 iG Cy / rot (9, 3) rot (9, 3) Ik d, tot (9. 3) grad hy dj 
® B B 


Cl, [ grad g,- 7 


rot (h, 4) 
+A [ga grad hy tot §)df + hy Ih d, J grad h, 
u 
Ber = J | 


und unter Benutzung der Relationen (3) für die Teilbereiche A os B läßt sich 
dafür schreiben: 


Le, (# -[ö > a gradh, - ds 


& 2 6G, we 
+4 fe > Cy (4 | z= k, d, | Ju grad h, - 
y 
=Clyl, baw. = [du gradh, - . ds 
(4b) 


y 


y 


2 
» | grad, - -d3 bzw. =—O° _ ) 
Y 


2 „ kennzeichnet die Differentiation von B aus gesehen nach außen. . 
Als eine Bestätigung für die Richtigkeit unseres Gleichungssystems wollen 


Energiebilanz für unser spezielles Problem zu formulieren. 
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ae: Der Energiestrom bei z = —oo muß gleich sein dem bei z = + uk Er ist bis 
auf einen rg Faktor gegeben durch den Realteil des Integrals: 

S, = J (Grr. x Of.) A af. 


Und die R <o} = R{S,> 0} lautet daher ausgeschrieben: 


— > 1,0, grad 


rot(h, 3) 
\; _ 


k, grad +5 tot (97 3) > 2 ke d; grad df | 


„grad g, l,c,gradg, + of 2 d,rot(h, 


rot(h, 3) 


O* — rot 5) 


En wobei 3 di = df gesetzt ist. Die Differentialoperatoren grad und rot sollen hier 
im folgenden durchweg als nur auf die Ebene senkrecht zur z-Achse bezogen 
verstanden werden. 


R Durch Ausmultiplizieren und Berücksichtigung der Orthogonalitäten: 
f grad 9% x gradh, df = 3/ (grad g, x gradh,)-3df=O 
G G 


J vot (9,8) tot (h, 3) df = 0. 
G 


— Man beweist sie, indem man die Integrale auf die Ausdrücke in (3) zurückführt — 
ergibt sich: 
[grad Ip X > tot 3) -df 
RK l,cygrad g, X 26 rot (9,3 )df 
ost tot(h, 3) x joe 


d, vot (hy 8) x I, ke de grad df 
* { G 
~ 


Erinnert man sich an die Beziehung: rot (h, 3) = gradh, X 3, formt die drei- 
fachen Vektorprodukte etwas um und wendet die Greenschen aa an, so stellt 
_ sich der Energiesatz in der endgültigen Form dar: 
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Er ist bis Die Summen erstrecken sich hier nur über die den homogenen Wellen ent- 

sprechenden Indizes, da die inhomogenen für den Realteil nichts liefern. 

; Wenn das Blendenproblem vollständig formuliert ist, müssen sich diese Rela- 

| tionen zwischen den Koeffizienten aus unserem Gleichungssystem folgern lassen. 
Um das zu zeigen, setzen wir die Existenz einer eindeutigen Lösung voraus und 
benutzen deren Eigenschaften in der Form: 


| — 2 c, rot (9,3) + 2 k,d, gradh, —( auf B 
Cl, grad g, -> l, cy grad gy + = d, tot (h, 3) 


| 
| | Wir begnügen uns mit der Behandlung bei einfallender Z-Welle, da sich der 
Fall einer H-Welle völlig analog erledigt. 


Trennung in Real- und Imaginärteil liefert: 


I cotG@,s) + |d,gradh, — > grad, = 0 (5a) 
df C y=] y=8+ 
2 auf B 
> rot 3) + |k,| dy grad hy + d, grad hh, = 0 (5b) 
len hier 6 y=] y=e +1 
bezogen r 
\1,| c, grad g, — c, grad 9, rot (hy 3) = grad gp (5c) 
=] 
auf A. 


r oo 
“ 2 ¢, grad 9, + 2 c, grad 9, tot (h,3) = C\l,| gradg, (5d) 


Die Zeichen v und ~ charakterisieren hier Real- und Imaginärteil der Koeffi- 
zienten. r bzw. s sind die Indizes der letzten homogenen Welle. 


führt — Wir integrieren das vektorielle Produkt der linken Seiten von (5a) und (5c) 
über @. Es ist dabei auch zulässig, von (5c) die rechte Seite zu nehmen. Denn im 
Bereich B, in dem die Gleichheit zwischen linker und rechter Seite nicht mehr be- 
> || hauptet ist, verschwindet ja der erste Faktor des Produktes wegen (5a). Benutzt 
(9,3) af man bei der Operation also einmal die linke und einmal die rechte Seite von (5c), 


S 


= 0. =r+1 (5e) 


e drei- 
) stellt Dieselbe Behandlung der Gin. (5b) a (5d) führt auf: we 
"> (Go. + — 2 20.06 
(5f) 


+28 (hy, hyo + > K3 (hy, = Lp Oy 


v=r+ 
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Abb. 2. Blendenform im 
rechteckigen Rohr 
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Durch Addition von (5e) und (5b) gewinnt man die oben abgeleitete Gestalt 


N bedeutet den Normierungsfaktor. 


hat damit gezeigt, daß er bereits im Ansatz enthalten ist, 
Auflösung des vorliegenden Gleichungssystems sind die uns 


bekannten Methoden ungeeignet. An einem speziellen Beispiel, das wir jetzt dis- 
kutieren wollen, werden wir aber die Eigenschaften der Matrixelemente genauer 
studieren und die Möglichkeit finden, eine inverse Form zu konstruieren, die das 
3 gesuchte Feld näherungsweise zu berechnen gestattet. 


II. Spezialisierung auf ein Rohr mit rechteckigem Querschnitt und schlitz- 
férmiger Blende und Ableitung des Projektionsmatrizenkalküls 


Der Querschnitt des Rohres werde als rechteckig mit den Seitenlängen a und} 
Das dem Problem angepaßte kartesische Koordinatensystem 
werde mit seinem Ursprung in den geometrischen Mittelpunkt der Blendenfläche 
gelegt, so daß = —+a/2 und y=-+b/2 die Ränder wiedergibt. Die Rohrwellen- 
funktionen lauten dann: P 


3) sın exilnm2 mit 


1 1 
3) cos mit = (6a) 


Die doppelte Indizierung empfiehlt sich wegen der Entartung der Eigenwerte: 


nen? m?n? a 
= 


- Lnm 2 
> 4 


Wir nehmen die Blende schlitzförmig an mit Schirmrändern parallel der y- 
Achse, gegeben durch 2 = + & (Abb. 2). Das Gebiet A entspricht demnach dem 


Koordinatenbereich § <|x | <a/2,—b/2 <y <+ b/2und 
B dem Koordinatenbereich |x| <&, —b/2<y< +6/2. 

Die y-Abhängigkeit der einfallenden Welle ist allen 
zur Beschreibung des Zusatzfeldes notwendigen Rohr- 
wellenfunktionen eigen, wie sich unmittelbar aus Ortho- 
gonalität und Randbedingungen schließen läßt. Dieser 
Umstand gewährt die Möglichkeit, auf die doppelte In- 
dizierung zu verzichten, wenn man vermerkt, daß der 
weggelassene Index für alle in Frage kommenden Eigen- 
funktionen derselbe ist. 

Aus (6a) folgert man weiter, daß für m = 0 nur 
H-Wellen existieren können. Die Blenden wandeln nicht 
um, d.h. bei einfallender H-Welle vom Index m = 
kann das durch Beugung und Reflexion entstehende 
Feld durch H-Wellen allein beschrieben werden. Darin 
liegt die eigentliche Spezialisierung unseres Problems. 


Aus z = — oo falle eine Welle vom Typus h,, ein. Das dazugehörige Eigen- 
system umfaßt alle h,, mit y= 1, 2,3... Die elektromagnetischen Feldkompo- 
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Weiter ist: / und K,= "2. 
Das Gleichungssystem (4b) vereinfacht sich fiir unser Beispiel wesentlich. Die 
durch Differentiation nach c, entstandenen Gleichungen fallen weg, und die c/ 
sind alle null. Wegen der Symmetrie des Schirmes verschwinden die Integrale 
über den Rand y, so daß schließlich bleibt: 


2 


n=1 n=1 


Setzen wir d, k, = y, und we daß über den Betrag von A noch frei ver- 


fügt und daher statt A auch ? -K geschrieben werden kann, so ergibt sich: 


a (aims Yn +2 (hin C (hms hi)s i 
oder: 
1 
a k, + Ak Bmn) Yn = 


n 

Die Matrix der &,„ heiße W, die der B,, B. Die x», sind Integrale des Produktes 
zweier normierter Orthogonalfunktionen eines vollständigen Systems, genommen 
über den Teilbereich A, die ß,„ über den Komplementärbereich B. Unser spezielles 
Funktionensystem ist vollständig im Sinne aller zu x = 0 symmetrischen und von 
y unabhängigen Funktionen. 

Die eben genannten Eigenschaften der Orthogonalität, der Normierung und 
der Vollständigkeit eines Funktionensystems ¢, bezüglich eines Gebietes @ be- 
liebiger Dimensionszahl und die Komplementarität der Bereiche A und B in dem 
Sinne, daß sie aus einem oder mehreren Teilgebieten bestehend zusammen das 
Grundgebiet G ausmachen, sind die einzigen Voraussetzungen, die wir der weiteren 
Betrachtung zugrunde legen. 

Wir wissen, daß jede im Lebesqueschen Sinn quadratintegrable Funktion 
Fourierkoeffizienten hat, deren Quadratsumme gegen das Integral des Quadrates 
der Funktion, genommen über das Grundgebiet, konvergiert. Umgekehrt ist nach 
Fischer-Riesz jedem Koeffizientensatz mit konvergenter Quadratsumme eine 
derartige Funktion äquivalent und zwar, abgesehen von einer Menge vom Maß 
null, in eindeutiger Weise. Die Frage, ob es zu einer gegebenen Funktion eine 


an 
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WERTEN nach den.g, gibt, durch die sie dargestellt wird, ist für die 

weitere Diskussion ohne Belang. Existiert das Lebesquesche Integral des Qua- 

drates einer Funktion, so wollen wir sie, wie gebräuchlich, eine Z?-Funktion nennen, 
Wir definieren die beiden Funktionen 


=, auf A =OaufA 
=0O auf B = Py auf B. 


Sie haben mit gy die Eigenschaft gemein, Z?-Funktionen zu sein. 
Ihre Entwicklungskoeffizienten sind: 


B 
= (Ps, Pu) _ = Box. 


Für zwei Z?-Funktionen f und g mit den Entwicklungskoeffizienten c, und d, 
gilt die Parsevalsche Gleichung: re 


= Sow. 

vl 
Demnach ist: 


co 
o=1 
bzw.: 
= (9; Yu)B = Pro Bou = Pru 


und ebenso: 


4) 
o=1 


i. 

Fr f 

A und ® sind demnach idempotente Matrizen, derart daß: 
B" = B AB = BA = 0. 


are 


Außerdem ist: 21+ 8 = € die Einheitsmatrix, denn us 


folgt: 
Kyu Bou = 

Wir wollen X und ® Projektionsmatrizen nennen. Ist nämlich € ein Vektor 
im Sinne Hilberts, dessen Komponenten c, die Entwicklungskoeffizienten einer 
L?-Funktion f sind, so ist der Vektor d, der aus € durch die Operation X entsteht, 
also d = UZ, einer Funktion f(4) zuzuordnen, die auf A mit f identisch und auf B 
identisch null ist. 

2 Zum Beweise stützen wir uns auf die Parsevalsche Gleichung und die Tat- 
sache, daß mit f auch f(4) eine L?-Funktion ist. Also: 


= d, = (7, 2°), = ou Cus 
u 


womit die Behauptung bewiesen ist. 
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Der Beweis fiir ® ist selbstverständlich derselbe, und die Projektionseigenschaft 


‘jst daher in bezug auf alle Z?-Funktionen gesichert. X und ® sind im Sinne Hil- 


berts beschränkte Matrizen, da sie jeden Hilbert-Vektor wieder in einen solchen 
überführen, oder weil die quadratische Form: 


für jeden Hilbert-Vektor € beschränkt ist. Era a 


Die Bilinearform EA®, wo € der L?-Funktion g entspricht, unterliegt einer 
ähnlichen Abschätzung und zwar auf Grund der Schwarzschen Ungleichung: 


Wir wenden uns voriibergehend wieder unserem speziellen Gleichungssystem 
m und schreiben es noch einmal in der Form an: 


(AK (7a) 


%, ist also jetzt ein Vektor im A-Raum, den die Operation X ungeändert läßt 
und die Operation ® zu null macht. Das heißt, Ba, =0. ist eine 
Diagonalmatrix mit den Koeffizienten 1/k,. Sie hat eine Reziproke $-! mit den 
Koeffizienten k,. Diese ist zwar nicht beschränkt. Aber das Produkt hat einen 
Sinn, wenn es auch, auf einen Hilbert-Vektor angewandt, eventuell aus dem 
Hilbert-Raum herausgeführt. Durch linksseitige Multiplikation von (7a) mit 
¥ folgt B K-18 Y= 0. Die quadratische Form 7 B 8-18 %Y ist aber dem Ausdruck 
Eok, äquivalent, wenn c, die Fourier-Koeffizienten der zu By gehörigen 
Funktionen sind. Daraus folgt: 8 ¥ = 0 wenn die k, alle gleiches Vorzeichen haben 
oder, wie in unserem Falle, alle das gleiche Zeichen abgesehen von der imaginären 
Einheit. Wir können daher statt (7a) auch schreiben: 


(AUR +A4B)y¥=—Cay. (7) 

Das spezielle Blendenproblem ist nun weit einfacher anzusetzen als auf dem 
Umweg über die Forderung des kleinsten Fehlerquadrats. Wir haben nämlich 
unsere Grenzbedingungen nur in der Sprache der Projektionsmatrizen zu formu- 
lieren. 

Dem magnetischen Transversalfeld seien die Entwicklungskoeffizienten y,, 
d.h. der Vektor y zugeordnet. Das elektrische Transversalfeld wird dann not- 
wendig durch & y gegeben sein, und die elektrische Komponente der einfallenden 
Welle mit bzw. C 6,. 


Die Grenzbedingungen lauten jetzt: wei 
ARF+CH)=O 


und: 


By =0. 


Multipliziert man die zweite Bedingung mit A und addiert sie zur ersten, so bekommt 
man die Gl. (7). 

Es sei nachträglich bemerkt, daß die Projektionseigenschaft der Matrizen auch 
für Funktionen mit komplexen Entwicklungskoeffizienten gilt. Die SR Ab- 
kitung wird dadurch kaum geändert. 
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Da sich für die spätere numerische Berechnung eine reelle Schreibweise als vor- 


____ teilhaft erweisen wird, trennen wir die Gl. (7) in Real- und Imaginärteil: 
oo oO, . 
- — + = % E Buv =9 
| k,| v=8+1 v=] re 
8 Y 00 Y ~ 
— — = * —— = Bur = 9. 


Wählen wir die einfallende Welle so, daß nur eine einzige homogene Rohrwelle 
möglich ist, also k®— K?>0 k®— K? <0 dann wird aus (8) nach einer kleinen 
Umstellung: 


~ 
00 Y §, y 


=O 
=> 


hay 
| 


eh) 
LS Ist das System (8b) eindeutig lösbar, so gilt dies auch für (8a) und umgekehrt. 
Daraus läßt sich das Verhältnis von Real- und Imaginärteil der Koeffizienten an- 
geben zu: 
—= ———— fiir alle». 


die noch zweckmäßig festgesetzt werden darf. Es sind daher die beiden Forderungen 


‘In der Amplitude C= C +7C steht eine komplexe Konstante zur Verfügung, 
zulässig: 


Sie bedingen im einzelnen die Relationen: 


~ > Yı Yı 
„=0 (=1-- ud C-CH=1. (9) 
Der Energiesatz — (C wA+C = fiihrt auf die Trivialitat: ik, 


wie wir erwarten müssen, da er schon in den Randbedingungen steckt. 
Die in erster Linie interessierende Größe, der Reflexionskoeffizient, ist als 
Funktion von %, bzw. d, darstellbar in der Form: 


= a, (10) 
yve+e V2yi—2yk +h 


Wie wir sehen werden, empfiehlt sich für die Berechnung die Einführung eines 


neuen Vektors € = — Dann wird aus (7): 
1 


; (UK + 2 ß, mit p= (Bu Biz, Bis - - -) 


6, = (1, 0, 0 ee >>. 


3 
|: 
iit 
g 
gan 
prü 
IV. 
7 vor 
| 4, 
der 
B- 
Un 
die 
in ¢ 
ni 
vol 
Ve 
80 
au 
i 
un 
na 
4 
17 
u 
ta 


(8) 


‚ohrwelle 
r kleinen 


(8 


(8b) 


gekehrt. 


nten an- 


E. Ruch: Einfluß einer Blende auf das Feld einer einfallenden elektromagnetischen Welle 263 


Wenn 1/A für A steht, führt die Multiplikation mit A.auf die Gleichung in der end- 
gültigen Gestalt: 


(11) 


K ist jetzt die Diagonalmatrix mit den Absolutbeträgen 1/|k,|. Der Uber- 
gang vom komplexen System (7) zum reellen System (11) ist, wie man leicht nach- 
prüft, nicht an die Existenz von nur einer einzigen homogenen Welle gebunden. 


IV. Die Auflösung eines Gleichungssystems mit unendlich vielen Unbekannten 
durch Konstruktion einer reziproken Matrix 


Nach diesen Zwischenbetrachtungen für das Blendenproblem soll wieder die 
vorher beschriebene Allgemeinheit des Orthogonalsystems gy, und der Bereiche 
A, Bund @ zugrunde gelegt sein. f und g seien vorgegebene L*-Funktionen mit 
den zugeordneten Hilbert-Vektoren € und d. f4) und g(®) bezeichnen die A- und 
B-Projektionen von f und g mit den Komponentensätzen @ = Uc und b = Bd. 
Unter D wollen wir eine Diagonalmatrix verstehen mit Elementen d,,, die alle 
positiv und von null verschieden sind und.unter einer festen Schranke K liegen. 

Wir stellen das Problem, eine Funktion X mit den Koeffizienten x, zu finden, 
die auf B mit g identisch und auf Grund einer vorgegebenen Funktionaloperation 2 
in eine Funktion 2 {X} übergeht, die auf A mit f identisch ist. Von der Operation 2 
müssen wir voraussetzen, daß sie sich als Diagonalmatrix mit den Eigenschaften 
von D darstellen läßt. 

Der Ansatz lautet: 


BE=b 


Versehen wir die zweite Gleichung mit dem Gewicht A und addieren sie zur ersten, _ 
so bekommen wir die Gleichung: ne 4 


aus der durch linksseitige Multiplikation mit A bzw. ® die beiden Ausgangsbe- 
dingungen entstehen. 
Wir schreiben auf Grund der Komplementarität A+B=€ stattdessen: 


und versuchen, eine reziproke Form dadurch zu gewinnen, daß wir ganz formal 
nach Potenzen von A—AUD entwickeln: 


— + (12a) 


Multiplizieren wir die einzelnen Glieder der Reihe aus. Das Produkt von X und 
WD ist nicht vertauschbar, es sei denn, es wird von rechts mit X multipliziert. 


Wir schreiben daher: 
(A —AAD) = (A —AAD) 


und diirfen jetzt in der ersten Klammer ausmultiplizieren, ohne uns um die Ver- 
tauschbarkeit zu kümmern, da jedes Produkt aus A und AD von rechts mit A 
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multipliziert ist. Daher ergibt sich: 


-(7 *) au) +... (AXDy. 


Beim Aufsummieren bekommt man nach dem Additionstheorem der Binomial- 


koeffizienten 
§ n 


die Partialsumme ©, in der Gestalt: 
J (12b) 
+{ +(7)4 
Der erste Klammerausdruck wird nach dem Binomialsatz: __ 
(ir = (€-AUD)". 
Der zweite ee die Umformung: 
1 
= - D) —€] D1 = — —(e-AUD)" 
so daß folgt: 
=|e- Dy" —= [(E-AUD)"—€] (12) 
Um die Existenz einer Grenzmatrix ‘ = lim ©, nachzuweisen, gehen wir zuerst 
n> 
von der Annahme aus, die d,, seien wesentlich von null verschieden, d.h. 
dy, >8>0 4 


mit einer festen positiven Zahl 4 für alle ». en 


Wählt man weiter A so, daß 2x <1, dann liegen alle Elemente 1 — A d,, zwi- 
schen 1—Ax>0 und 1—/Aöd<1 und die Operation €—AD, angewandt auf 
irgendeinen Hilbert-Vektor x verursacht eine Verminderung seines Betrages 
mindestens um den Faktor 1 —/ 6. 

Da die Matrix X ebenfalls den Betrag eines Vektors verkleinert oder höchstens 
gleich läßt, gilt die Ungleichung: 


|z(M—AMD)" z| = |< (1-Aö). 


Damit ist eine Majorante für die Reihe der Absolutbeträge der quadratischen 
Formen gefunden, die für alle Vektoren x des Hilbert-Raumes gleichmäßig kon- 
vergiert. Die Existenz einer beschränkten Grenzmatrix y ist also erwiesen. 

Bezeichnen die d,, eine Nullfolge, so ist diese Schlußweise nicht mehr möglich. 
Wir können aber dann unsere Matrix D in eine Matrix D’ abändern in der Weise, 


fr a einem Index » = N ab alle nachfolgenden Elemente dy,» durch dy et 
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setzt werden. Diese Abänderung ist beliebig klein zu machen, da N nur endlich 
sein muß, im übrigen aber beliebig groß sein kann. Mit der Matrix D’ gilt die 
eben vorgeführte Abschätzung. Daß der Grenzübergang © = lim G, für Null- 
n>0o 
folgen auch ohne jenen Umweg sinnvoll bleibt, liegt nahe. Wir werden uns 
damit begnügen, dies an einem leicht zu behandelnden Spezialfall zu bestätigen. 
Die weitere Diskussion setzt die Existenz von © voraus. Um zu beweisen, daß 
der Ausdruck: 
lim \e- AY D)" — (13) 


n> 
tatsächlich die inverse Matrix zu ® + AY D ist, zeigen wir, daß Links- wie Rechts-, 
multiplikation von y mit 8 + AA D auf die Einheitsmatrix € führt, was nach einem 
Satz der Gruppentheorie die Eindeutigkeit des inversen Elementes garantiert. 
Dafür ist eine kleine Vorbemerkung nötig: Die Existenz einer beschränkten Grenz- 
matrix lim (Ee—ANU D)" läßt die adäquate Formulierung zu: Die Vektorfolge 


n>00 
= (€ -AQND)rZ strebt gegen einen Grenzvektor Zo für alle Hilbert-Vek- 
toren x. Dazu ist notwendig, daß 


> > 


oder x, —(€ —AUD) x,—>0. Daraus folgt aber: Dax, — 0. 


n>00 


ta) dr 
lim (€—A%D)"—€] D> 


nur den A-Raum auf sich transformiert, also links- wie rechtsseitige PEN 
mit ® die Nullmatrix ergibt. Unter Verwendung dieser beiden Eigenschaften folgt: 
(B+ayD) —AUD)" — [(E-AA — 


denn B(€E— AU D)" = B gilt unabhängig von n. 3 
Die Linksmultiplikation ergibt: 


n>0 


= lim (B+AUAD) — [(E-AAD)"—E] 
n>00 
= lim 


n>00 


= lim 


n>00 


LD)? 19 4 
inomial- 
(12b) 
> 
(12) Weiter sieht man aus der Darstellung (12b),daß der zweiteTerm von S,nämlich: 
. 
ir zuerst 
Br 
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In dem ersten Produkt (€ —AYAD)"AUAD ist AD mit der Potenz vertausch- 
bar, im zweiten (E— AU D)" A, ist die Vertauschbarkeit zu erreichen, indem man 


von rechts mit D D- multipliziert. Dann wird: ‘i 
lim ADD =AD lim (€—-AAD)" D4. 
n-> 0O n>00 


Diese beiden Ausdriicke sind aber nach der Vorbemerkung null. Also ergibt sich 
auch hier, wie bei der Linksmultiplikation die Einheit €, und es ist gezeigt, daß & 
im Hilbert-Raum eine eindeutige reziproke Matrix zu 8 + AY D ist. 

Fiir Nullfolgen ist eine kleine Abweichung zu erwarten. Wir werden das an 
einem Beispiel sehen. 

Nehmen wir an, A sei das ganze Grundgebiet G. Dann vereinfacht sich die 
Gleichung 


ADz=ia 


Die reziproke Matrix heißt also: + D>. 


: Damsibe muß natürlich auch mit unserer Matrix © herauskommen. Wir sehen: 


(l—Ady)" 0 0 
lim (€—AD)" = lim 
n>00 0 0 (1 


und aus © wird +D- 1 wie oben. 


Dabei ist es RR zulässig, daß sie d,, eine Nullfolge darstellen. Wir können 
nur nicht mehr behaupten, daß © eine beschränkte Matrix ist, sondern müssen 
formulieren: sie ist sinnvoll im Sinne einer quadratintegrablen Lösung, wenn sich 
die Operation D-! auf die rechte Seite anwenden läßt, ohne aus dem Hilbert- 
Raum herauszuführen. 

In der allgemeinen Form von © tritt die nicht beschränkte Matrix D-! nur beim 
zweiten Term auf, der, wie wir schon festgestellt, den A-Raum allein auf sich 
transformiert. Die Einschränkung, die der rechten Seite aufzuerlegen ist, um 
eine quadratintegrable Lösung erwarten zu können, ist wesentlich komplizierter 
und soll hier nicht untersucht werden. Die Tatsache, daß bei Abänderung der 
Matrix D in D’ die Schwierigkeiten entfallen, enthebt uns der Notwendigkeit 
einer Untersuchung, wenn die Existenz einer Lösung aus physikalischen Gründen 
feststeht. 

Nach diesen abstrakten Entwicklungen können wir unser spezielles Gleichungs- 
system (11) für die Blenden im rechteckigen Rohr bequem diskutieren und auf- 


lösen. 
ir schreiben es noch einmal an: & 


% z charakterisiert das magnetische Transversalfeld auf B und AK x das elektri- 
auf 
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Da die rechte Seite nur auf B von null verschieden ist, lautet die reziproke 


Matrix X 


2 


V. Numerische Berechnung des behandelten Spezialfalles 


Zur numerischen Berechnung nehmen wir eine endliche Anzahl von Gleichungen 
mit ebenso vielen Unbekannten. Wir wollen uns mit 16 begnügen. Wir wenden die 
Operation €—AMK auf den Vektor —ß, an und bezeichnen das Ergebnis als 
erste Näherung x). Dieselbe Operation, auf den Vektor 2“ angewandt, gibt uns 
dann die Näherung z® usw. Da die vorangegangenen Überlegungen bezüglich 
der Projektionseigenschaften der Matrizen bei der Beschränkung auf endlich viele 
Elemente nicht mehr genau stimmen, wird sich bei jeder Operation ein Fehler ein- 
schleichen, der die Häufigkeit ihrer sinnvollen Anwendung beschränkt. Das Ver- 
fahren konvergiert also nicht streng gegen die gesuchte Grenzfunktion. Aber 
wenn wir uns vergewissern, daß der durch die Endlichkeit bedingte Fehler jeweils 
klein ist, können wir eine bestimmte Anzahl von Operationen durchführen und 
erwarten, daß das so gewonnene Ergebnis mit befriedigender Genauigkeit die Lö- 
sungsfunktion darstellt. 


Die schlitzförmige Öffnung mache die Hälfte des Gesamtquerschnitts aus. Die 
Schirmränder liegen also bei &=+.a/4. Dann sind die 


a 
4 


2 
1 Bon ==> sin va (= +5)sinua (+ 


bol m 


Mit der Substitution z (= — 5) =z bekommen wir: 


Die elemenatre Integration führt fiir vy + u auf: Fin a eo kr 


3n 


- 4 
Einsetzen der Grenzen unter Benutzung der trigonometrischen Beziehung : 


sin 3a — sina = 2 cos 2a 
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Wir tien, 4 


für v + u ungerade wird das Integral null 


fiir ungerade erhält man—(—1)* ?. 
Wir können daher schreiben: oh 


1 
(—1) * 2 (—1) * 2 
Wenn Erz für einen Index, der ungerade oder durch vier teilbar ist, die Zahl null 


und für einen Index, der gerade aber nicht durch vier teilbar ist, die Zahl eins be- 
deutet. Die Auswertung für » = u ist wesentlich einfacher und ergibt: 


+1 


1 für und gleich null für vy sei. Die Matrix 


Der Vektor — B, auf der rechten Seite hat die Komponenten x en ao 


= 7 (Bu; Bis 


Er hat, wie man sieht, die Eigenschaft, daß die geraden Komponenten jeweils 
null sind. Die Gestalt der Matrix 2{ und damit auch B und A X ist derart, daß die 
durch die Operation €—AUK hervorgehende nächste Näherung diese Eigen- 
schaft bewahrt. Dies hat die erfreuliche Folge, daß wir mit weniger Rechenarbeit 
auskommen. Denn wir können jetzt für unsere Matrizen und Vektoren auf alle 
geradzahligen Komponenten verzichten. 
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‚fiehlt es sich, die Matrix X K zu bilden und die Operation 
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Die ursprünglich auf sechzehn Zeilen reduzierten Matrizen vereinfachen sich 
damit praktisch zu achtzeiligen. Wir setzen daher in Zukunft für 2»—1 den 
Index v. Eine triviale Rechnung zeigt uns, inwieweit die Projektionseigenschaft 
der Matrix 2{-noch zutrifft. Kombinieren wir die &,, mit den &,,, so müßte a2, 
also 1/ herauskommen, wenn alle Elemente, d.h. unendlich viele, berücksichtigt 
wären. Bei Beschränkung auf die ersten acht Komponenten ergibt die Rechnung 
e- , also einen Fehler von etwa 3,7%. Für größere Indizes, beispielsweise agg ist 


der Fehler natürlich erheblich größer. 


Die Frequenz der einfallenden Welle wählen wir so, daß k= En wird. Dies 
entspricht einem Verhältnis = = = , wenn mit w, die Grenzfrequenz bezeichnet 


ist, unterhalb der ıeelle Wellen nicht mehr existieren. Für die Eigenwellenzahlen 
gilt dann mit der neuen Indizierung 


k= 


und es zeigt sich, daß k, reell, k, bereits imaginär ist, d. h. nur die homogene Welle 
vom niedrigsten Typ existenzfähig ist. Bei der numerischen Behandlung emp- 


vorzunehmen. Der Vektor AK x") gestattet dann zu beurteilen, inwieweit die 


Forderung des Verschwindens von K z im A-Raum erfüllt ist. Das Verfahren hört 
damit auf, daß der sich infolge der Endlichkeit der Matrizen einstellende Fehler 
die B-Projektion des Ausgangsvektors wesentlich zu verändern beginnt. 


Damit der Prozeß für unendlich viele Glieder streng konvergiert, muß gelten: 


Innerhalb dieser Grenzen ist das A willkürlich wählbar. Mit kleinem A Zr der 


Fehler kleiner, dafür werden aber mehr Näherungen nötig sein. Im Falle — ~ Te] * wl 


wird man mit weniger Rechenarbeit denselben Effekt erzielen, weshalb wir uns. 
für die letztere Wahl entschließen und A = az setzen, denn: 


Min |p | = || = 3; 2,286. 


Mit — A auviuien es sinnvoll, wie aus den Abb. 3 und 4 hervorgeht, sich 
mit der vierten Näherung zu bescheiden. Wir geben die BE; zweite und vierte 


| | 4 
' 
— | 
5 
0 
| 
3 
n jeweils 
, daB die Br 
>» Eigen- 
enarbeit 
auf alle 


. 1950 


at xf) 
—2,571 | —2,122 —1,948 | 2,236 —1,151 —0,949 —0,871 
+1,000 +1,723 +1,966 5,196 +0,193 +0,332 +0,378 
+0,333 +0,481 +0,421 9,539 +0,035 +0,050 +0,044 
—0,333 —0,694 | —0,906 | 13,675 —0,024 —0,051 —0,066 
—0,200 —0,336 —0,389 | 17,748 —0,011 —0,019 —0,022 
+0,200 +0,429 +0,581 | 21,795 +0,009 +0,020 +0,027 
+0,147 +0,254 +0,287 | 25,826 +0,006 +0,010 +0,011 
—0,147 —0,303 —0,406 | 29.850 ‘5 —0,005 —0,010 —0,014 
Für die zeichnerische Darstellung führen wir wieder den Vektor y= (2 + 4) || 


ein. ¥ gibt uns das zusätzliche magnetische er und § % das elektri- 


sche, wenn die einfallende Welle die Amplitude C = 72 + i(a- ir ie) besitzt. 


Yı Ye 


Die erste Komponente von y= (2 = ‘) ist reell, die übrigen sind imaginär. 


Esempfiehlt sich, zur graphischen Veranschaulichung statt RY Ky y= Gz El er) 
zu nehmen. Er gibt einfach die skalare Differenz von Imaginär- und Realteil des 
magnetischen Zusatzfeldes und ist daher auf die Amplitude C’ = C-(=-=4 
zu beziehen. Durch diese Überlegung ist noch einmal der Unterschied der beiden 
Gin. (7) und (8b), also zwischen (A & + 2B) ¥ =— Ca, und (AK+AB)Y 


erläutert. 


Der Vektor K Y stellt genau so wie & Y das elektrische Feld dar, nur zu einem 
anderen 


= 


Abb. 3. Magnetisches Transversalfeld. Die Kurven 0, 2, 4 bedeuten der Reihe nach 
die nullte, zweite, vierte Näherung 


Die Abb. 3 und 4 geben das magnetische Feld Y und das elektrische K Y. Die 
fein ausgezeichnete Sinuslinie 1 in Abb. 4 ist die negative Amplitude plus eins. 
deren stark ausgezogener Teil zusammen mit der x-Achse im Intervall von —a/4 
bis +a/4 in Abb. 3 die Forderung veranschaulicht, die an K Y bzw. % durch die 
Randbedingungen gestellt sind. Die punktierte Kurve 0 ist die nullte Näherung, 
die bei Berücksichtigung von unendlich vielen Fourierkoeffizienten in Abb. 4 
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auf A mit der abgeschnittenen Sinuslinie und auf B mit der x-Achse zusantmen- 
fallen müßte. Man sieht aus der Abbildung, was man von der Lösung bei Be- 
schränkung auf 16 Fourierkoeffizienten an Genauigkeit zu erwarten hat. Die ge- 
strichelten Kurven 2 und die ausgezogenen 4 geben die zweite und vierte Näherung. 
Man erkennt an Abb. 3 das Anwachsen des magnetischen Feldes auf dem Schirm 


te 


Abb. 4. Elektrisches Transversalfeld. Die Kurven 0, 2, 4 bedeuten der Reihe nach 

die nullte, zweite, vierte Näherung. Kurve 1 ist die einfallende Welle mit entgegen- 
gesetztem Vorzeichen 


gegen den Rand zu, das, wie wir wissen, am Randy gegen unendlich gehen muß. 
Die geforderte Übereinstimmung zwischen negativer Amplitude und Zusatzfeld 
des elektrischen Vektors ist mit der für 16 Fourierkoeffizienten möglichen Genauig- 
keit durch die vierte Näherung zufriedenstellend wiedergegeben, wie aus Abb.4 9g 
hervorgeht. Die starke Oszillation des magnetischen Feldes in der Öffnung ist 
durch die schlechtere Konvergenz der Lösung gegenüber der Ausgangsfunktion 
verständlich. Außerdem haben wir ja schon bei der Beurteilung der Projektions- 
eigenschaften der Matrizen hervorgehoben, daß die höheren Entwicklungskoeffi- 
zienten ungenauer zu erwarten sind als die ersten. Der erste Entwicklungskoeffi- 
zient aber ist es gerade, der die wesentlichen Verhältnisse, nämlich Reflexions- und 
Durchlässigkeitsvermögen der Blende, wiedergibt. Man errechnet sie nach der 
Formel (13) zu: 
= 0,68. 


Der reflektierte Energiebetrag ist also: R = (0,68)? = 0,46 ae, a1 
und die durchgelassene Energie: D= 0,54. ie na 


VI. Diskussion des allgemeinen Blendenproblems 


Wir wollen noch einmal auf das allgemeine Blendenproblem, also auf die 
Gln. (7a) und (7b) zurückkommen. Es unterscheidet sich im wesentlichen dadurch 
von dem behandelten Fall, daß E- und H-Wellen gleichzeitig auftreten. Aus 
den Relationen (3) folgt, daß das vollständige System der Funktionen grad 9, 
und rot (A, 3) orthogonal ist. Normiert man dieses Funktionensystem und ordnet 
die Gleichungen und in den Gleichungen die einzelnen Glieder nach steigenden 4 
Eigenw erten der E- und H-Wellen gleichzeitig, dann läßt sich unser System (4a) 
im Projektionsmatrizenkalkiil genau so schreiben wie das spezielle System (7). — 
Man weiß, daß zwischen zwei Eigenwerten der schwingenden Membran mit ein- 
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gespannten Rändern ein Eigenwert der schwingenden Membran mit freien Rand 
liegt und umgekehrt. — Das allgemeine Problem ist demnach genau so behande 
wie das spezielle. Die Ähnlichkeit der durch Differentiation nach den c¥ di 
seits und nach den d* andererseits gewonnenen Gleichungen legt die Vermutill 
nahe, daß wenigstens in Spezialfällen durch geeignete Umformung zwei Syst¢ 
aufgestellt werden können, die nicht miteinander gekoppelt sind und von dem 
die Lösung des einen genügt, um die Lösung des anderen anzugeben. Diese Möglig 
keit würde implizit eine Aussage über das Umwandlungsverhältnis enthalten ı 
eine erhebliche Vereinfachung der Rechenarbeit liefern. 

Die Allgemeinheit der mathematischen Fragestellung und die Möglichkeit i 
Beantwortung nach einer Methode, die ein erträgliches Maß an Rechenarbeit ve 
langt, gibt Veranlassung, ähnliche Probleme in der Sprache des Matrizenkalkä 
zu formulieren. Es liegt im Wesen des Ansatzes und der Methode, daß dabeiä 
zweiter Linie die Verhältnisse an den Grenzflächen der Teilbereiche wiedergegebe 
werden — sie werden durch die hohen Fourierkoeffizienten beschrieben — währe 
in erster Linie Größen wie Reflexionsvermögen und Durchlässigkeit beim Roi 
wellenproblem, also allgemein solche, die wesentlich von den ersten Fourig 
koeffizienten abhängen, bestimmbar sind. ; 


Herrn Prof. Dr. G. Hettner, an dessen Institut die vorliegende Arbeit angg 
fertigt wurde, möchte ich für freundliche Unterstützung an dieser Stelle meine 
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der Arbeit und Herr Dr. F. Wenzel war so liebenswürdig, den Teil über dm 
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